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Laser-Induced  Line  Narrowing  h'ifeci.  in  Coupled 
Doppler-Broadened  Transitions 
M.  S.  I 'eld  and  A.  Javan 


ABSTRACT 


The  line  shape  of  n  Doppler-broadened  transition  is  dramatically 
altered  by  the  presence  of  a  laser  field  resonating  with  a  second  Doppler- 
broadened  transition  sharing  a  common  level:  two  narrow  resonances  of 
different  widths  appear  superimposed  upon  the  oroud  background  signal 
at  frequencies  symmetrically  located  about  the  corresponding  line  center. 
The  effect  has  already  found  application  in  a  number  of  seemingly  differ¬ 
ent  though  intimately  related  studies,  including  high  resolution  h.  f.  s. 
and  isotope  shift  determinations.  The  theory  of  the  effect  is  developed 
with  reference  to  these  applications.  The  treatment  is  formulated  in 
terms  of  transition  rates  induced  by  two  classical  fields  resonantly  inter¬ 
acting  with  a  pair  of  coupled  Doppler-broadened  transitions  of  arbitrary 
frequencies.  The  perturbation  approach  adopted  is  valid  for  one  field 
fully  saturating  its  transition;  the  resulting  line  shape  expression  exhibits 
important  power  broadening  effects.  This  approach  is  equivalent  to  the 
familiar  density  matrix  formulation,  which  is  also  presented.  Various 
features  of  the  resulting  expression  are  discussed  in  detail  as  they  apply 
to  two  precision  spectroscopic  applications,  Mode  Crossing  and  Sponta¬ 
neous  Emission  Line  Narrowing.  The  connection  with  previous  work  is 


also  discussed. 


FIGURE  1:  Spectrum  Of  Lino  Narrowing  Induced  In  The  Doppler  C> 
(Or  Attenuation)  Profiles  Of  (a)  Laser  And  (b)  Coupled  Transitions 


I.  INTRODUCTION 


It  is  '.veil  know ^  that  the  overall  gain  profile  of  a  Doppler-broadened 
laser  transition  is  dramatically  influenced  by  the  presence  of  the  laser  field. 
This  may  be  demonstrated  by  scanning  the  gain  profile  with  a  weak,  mono¬ 
chromatic  probe  field  eolinear  with  the  laser  field;  As  the  probe  field  is 
tuned  through  the  transition,  two  identical  sharp  decreases  in  gain  appear 
superimposed  upon  the  broad  lineshape,  one  at  the  laser  frequent y  and  one 
symmetrically  located  on  the  opposite  side  of  the  atomic  line  center. 

[Fig.  1(a)  ).  These  resonant  decreases  occur  because  the  standing-wave  field 
within  the  laser  cavity  selectively  interacts  with  atoms  whose;  velocities 
Doppler- shift  one  of  its  travelling-wave  components  into  resonance.  This 
produces  changes  in  the  laser  level  populations--  an  increase  in  the  lower 
level  population  and  a  decrease'  in  the  upper  level  population--  over  two 
intervals  symmetrically  located  about  the  renter  of  the  velocity  distribution. 
These  changes  reflect  themselves  in  depletions  in  the  gain  profile-  over  the- 
corresponding  frccjue-nery  inte  rvals.  The-  exte  nt  of  these  intervals  is 
dctrrmine-el  bv  the  natural  widths  (or,  more-  ge-nc-rnlly,  the*  homoge-ne-ous 
widths)  of  the  atomic-  transitions. 

In  the*  fore-going  discussion  it  is  assumed  that  the  laser  field  is  de¬ 
tuned  from  w9,  the  atomic  center  frequency.  As  its  frequency  approaches 
within  a  natural  width  of  u9,  the  two  resonant  decreases  merge  into  a  single 
one*.  However,  let  us  primarily  consider  those  cases  in  which  the  laser  fre¬ 
quency  is  detuned  and  the  change  signals  are  well  resolved. 
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The  line  shape  of  a  second  Doppler-broadened  transition  formed  by 
either  of  the  levels  of  the  laser  transition  and  a  thn  tl  level  is  also  consider¬ 
ably  altered  by  laser  oscillation:  Scanning  as  before  the  gain  for  attenuation) 
profile  of  this  transition  with  a  weak  probe*  field,  one  again  finds  two  sharp 
resonances  at  frequencies  symmetrically  located  about  the  corresponding 
line  center.  These  change  signals  differ  from  those  of  the  laser  transition 
in  one  remarkable  aspect:  one  can  be*  considerably  narrower  than  the  other. 
[Fig.  1  (b)|  .  Wo  shall  refer  to  the  widths  of  the  broad  and  narrow  change 
signals  as  T,,  and  rv,  respectively .  For  example,  if  the 
center  frequency  of  the  coupled  transition  (0  -  1 )  is  close  to  that  of  the  laser 
transition  (0-2)  then,  for  a  weakly- saturating  laser  field,  r^.  =  Yj  +  y2  anfl 
FR  =  Yj  +  Y2  +  ^  Yq.  Yj  drra>r  rato  l°vxd  j-  In  comparison,  the 

change  signals  of  t lio  laser  transition  are  each  of  w'idth  Tj  =  Tq  +  Y  2‘ 

Similar  linoshape  features  would  be  observed  in  the  closely  related 
situation  in  which  the  fluorescence  from  either  of  the  laser  levels  to  a  third 
level  is  monitored  along  tilt-  laser  axis.  Note  that  the  resulting  spontaneous 
emission  socctrum  directly  follows  the  spectrum  of  emission  stimulated  by 

a  weak  probe  field  tuned  through  the  coupled  transition  when  the  lower  level 
population  of  the  coupled  transition  is  ignored.  As  an  illustration,  suppose  the 
laser  field  were  tuned  to  the  low  frequency  side  of  its  Doppler  profile:  For  emission 

originating  in  its  upper  level,  the  laser-induced  change  signals  would  appear 
as  resonant  decreases  of  widths  T^  below  and  1^  above  with  the 


utllw.. 
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fluorescence  center  frequency.  In  contrast,  emission  originating  in  the 
lower  laser  level  would  result  in  resonant  increases  with  the  positions 
of  the  broad  and  narrow  change  signals  interchanged.  A  further  note¬ 
worthy  distinction  is  the  differing  radiative  origins  of  the  change  signals: 
in  ’he  former  case  they  result  primarily  from  double- quantum  transitions, 
while  in  the  latter  case  they  are  primarily  due  to  single-quantum  tran¬ 
sitions.  This  important  distinction  is  elaborated  below. 

Since  r^,  and  are  generally  much  narrower  than  the  Doppler 
widths  we  shall  refer  to  this  effect  as  "Laser- Induced  Doppler  Line 
Narrowing". 

As  has  been  noted,  the  stimulated  and  spontaneous  versions  of 

Laser- Induced  Doppler  l  ine  Narrowing  are  different  manifestations  of 

the  same  basic  quantum-mechanical  effect.  Nevertheless,  important 

applications  of  these  versions  bear  little  res  mblnnee  to  one  another  on 

the  surface.  A  major  purpose  of  this  paper  is  to  relate  the  basic  effect 

to  those  intimately  connected  though  seemingly  different  applications. 

In  some  of  these  L,  »  L.,  causing  the  differing  width  characteristic.' 
ij  IS 

to  reveal  themselves  in  particularly  striking  ways. 

Several  puibieations  dealing  with  details  of  the  line  shape  in  various 

(*>  _  o  \  (2) 

special  eases  have  appeared  previously.  '  The  initial  presentation  by 

Sehlos-bf  r^  and  .lav  n,  a  quantum-mechanical  analysis  of  the  third-order 
polarization  induced  by  two  classical  fields,  demonstrated  the  applicability 
of  the  narrow  resonance  to  high-resolution  studies  of  closely-spaced 
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Doppler-broadened  laser  transitions  (l'ig.  2(b)).  The  latter  treatment  also 

spplies  to  a  cascade  system  |l;ig.  2(c))  in  which  the  middle  level  lies  about 

halfway  between  the  upper  and  lower  levels.  Subsequently,  Nothin,  Hautian 

(4) 

and  Keoktistov  presented  a  nuanti/.ed  field  eairulation  whic  h  desc  ribed 

the*  spontaneous  e*mission  speed  rum  arising  from  one  of  the*  levels  of  a 

(ft) 

weakly-saturated  laser  transition.  A  recent  discussion  by  Holt,  formu¬ 
lated  on  the  basis  of  two-photon  transitions  induced  by  the  laser  field, 
analyzed  the?  frequency  profile  of  the-  spontaneous  emission  arising  from  the 
lower  laser  level  (Fig.  2(c)),  This  treatment  neglects  single-quantum 
events  and  would  not  apply,  for  example,  to  the  ease  of  spontaneous  emis¬ 
sion  arising  from  the  upper  laser  level.  The  treatment  of  Ref.  4,  how¬ 
ever.  is  valid  in  either  ease.  A  subseepuont  letter  ^  by  the  present  authors 
analyzed  the  spontaneous  and  stimulated  versions  of  Laser-Induced  Doppler 
lane  Narrowing  experiments  by  means  of  a  classical  field  approach, 
emphasizing  their  close  relationship;  the  theoretical  treatment  took  into 
account  the  influence*  of  both  single*-quantum  and  double-quantum  transitions, 

(<j) 

and  included,  in  addition,  intensity-dependent  line;  broadening  c*ffeets. 

(7' 

The*  present  pape  r  is,  in  par  t,  an  elaboration  of  that  letter,  *  and 
contains  additional  detailed  discussions.  The  treatment  is  formulated  in 
tea  ms  of  transition  rates  induced  by  two  classical  fields  resonantly  inter¬ 
acting  with  a  pair  of  coupled  Doppler-broadened  transitions  of  arbitrary 
frequencies.  The*  method  of  analysis  is  an  extension  of  the  one'  u,adopted 
some  time  ago  in  ealc-ulating  the  line*  shape  details  of  u  three  level  maser. 


In  this  approach  two  distinct  processes  emerge:  the  first  a  double-quantum 
tr  utsition,  involves  the  exchange  of  a  photon  with  each  of  the  two  applied 
fields:  the  second,  an  inherently  single -quantum  act.  includes  the  influence 
of  one  field  on  the'  rate  at  which  single-quantum  transitions  are  induced  by 
the  other  field.  This  distinction  is  not  apparent  in  the  usual  density  -  mat  fix 
formalism  in  whieh  the  induced  pole ri '.ation  is  calculated.  The  two  ap¬ 
proaches  are,  of  course,  equivalent:  their  connection  will  be  cl  irified 
below.  Moreover,  the  theoretical  approach  adopted  is  not  restricted  to  a 
third-order  polarization  calculation,  and  is  valid  for  one  field  fully 
saturating  its  transition,  the  resulting  lino  shape  expression  exhibits  im¬ 
portant  power  broadening  effects.  We  ire  able  to  obtain  such  an  expression 
because  in  the  applications  discussed  ..ere  it  is  generally  sufficient  to  con  ¬ 
sider  the  standing- wave  laser  field  to  be  detuned  from  the  center  of  its 
Doppler-broadened  gain  profile.  Then  its  travelling-wave  components  do 
not  couple  to  each  other  and.  consequently,  may  be  treated  independently. 
Furthermore,  major  applications  and  their  important  physical  features  are 
adequately  described  by  considering  the  uc  ak  field  gain  of  one  transition  as 
influenced  by  a  fully-saturated  coupled  transition;  the  approach  may  be 

readily  extended,  if  necessary,  to  include  higher  order  effects. 

(4 ) 

The  approach  of  Nothin,  Hautian  and  Feoktistov,  which  calculates 
the  laser- induced  spontaneous  emission  spectrum,  is  complementary  to  ours. 
That  treatment  is  formulated  considering  both  radiation  fields  in  quantized 
form.  The  actual  calculation  is  somewhat  simplified,  however,  by  con¬ 
sidering  the  laser  field  in  its  classical  form,  a  procedure,  of  course, 


justifiable  for  radiation  oscillators  in  states  of  high  excitation.  The  spon¬ 
taneous  emission  field,  however,  is  kept  in  its  quantized  form.  The 
introduction  of  field  quanta  leads  naturally  to  a  mathematical  development 
considerably  different  from  our  classical  field  approach.  The  ('(imputa¬ 
tions  are  rather  involved,  and  relevant  results  arc  obtained  only  for  a 
\v  akly-saturating  laser  field.  As  discussed  below  the  latter  results 
agree  with  the  weak-saturat ion  limit  of  those  obtained  here.  In  the  dis¬ 
cussions  cf  Ref.  4  the  connection  with  experimental  observations^*’  *^  ^ 

(2) 

and  theory  of  an  earlier  stimulated  version  of  the  effect  has  not  been 
made.  This  connection,  which  is  not  apparent  on  the  surface,  is  clarified 
below. 

The  calculation  of  the  Doppler-br  oadon  >d  response  is  carried  out 
in  i'art  II.  Section  IIA  considers  the  interaction  of  two  monochromatic 
classical  travelling- wave  fields  v.  ill)  a  group  of  three  level  atoms  moving 
with  fixed  velocity.  One  field  is  assumed  to  fully  saturate  the  transition 
with  which  it  resonates;  the  second  field,  assumed  weak,  probes  the  coupled 
transition.  The  emitted  power  at  the  probe  frequency  is  calculated  ir.  terms 
of  transition  rates,  as  indicated  above.  The  subsequent  average  over  atomic 
velocities,  which  leads  to  the  Doppler-broadened  travelling-wave  response, 
is  outlined  in  Section  IIB,  together  with  a  discussion  of  the  resulting  line 
shape.  The  extension  of  the  travelling- wave  analysis  to  standing- wa  /e 
applications  is  presented  in  Section  IK’,  which  also  discusses  the  specifically 


standing- wave  c.fects  arising  when  the  intense  field  is  tuned  to  its  atomic 
center  frequency.  Wherever  possible,  the  presentation  of  detailed  alge¬ 
braic  m.'.iipulations  and  proofs  has  been  deferred  to  Appendices. 

Part  III  discusses,  by  example,  important  features  cf  Laser- 
induced  Doppler  Line  Narrowing.  Section  I1IA  examines  the  frequency 
dependence  of  the  atomic  response  prior  to  Doppler  averaging.  Section  IIIB 
continues  by  discussing  the  extension  of  the  line  narrowing  effect  to  the 
spontaneous-emission  version  mentioned  above,  explaining  the  connection 
with  earlier  formulations.  In  Sections  IIIC  and  IIID  two  important  appli¬ 
cations  are  examined  in  detail:  The  first, a  technique  involving  two 
classical  fields,  enables  structure  of  Doppler-broadened  systems  with 
closely  spaced  levels  to  be  measured  with  great  accuracy.  This  tech¬ 
nique  has  already  been  employed  in  measuring  hvperfine  structure^12)  and 

paramagnetic  properties'2^  of  several  excited  atomic  levels.  The 
(14) 

second,  based  upon  the  spontaneous  emission  version,  has  been  utilized 
in  isotope  shift^'^  and  linewidth  paramete  measurements  in  No. 

The  results  of  this  paper  are  also  directly  applicable  to  the  extra¬ 
ordinary  behavior  of  atomic  Oxygen  fine- structure  laser  oscillations  at 
8446$,  already  briefly  discussed  in  Kefs.  8,  7  and  ID.  In  that  case 
other  physical  processes,  entirety  unrelated  to  the  present  discussions, 

arc  of  great  importance.  We  prefer  to  discuss  these  together  in  a  separate 
(90) 

publication  '  which  will  utilize  expressions  derived  below. 
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II.  DOPPLER -BROADENED  GAIN 

\\  f*  now  proceed  to  calculate  th<  interaction  of  a  weak,  mono¬ 
chromatic  probe  field  with  one  of  the  transitions  of  a  Doppler-broadened 
three  level  system,  as  influenced  by  a  saturating  field  resonating  with 
the  coupled  transition.  In  order  to  describe  applications  within  a  Fabry- 
Perot  cavity  one  must  consider  the  possibility  of  fields  in  the  form  of 
standing  waves,  as  well  as  travelling  waves.  For  reasons  given  below, 
the  major  features  of  important  standing- wa ve  applications  may  be 
understood  by  analyzing  the  ease  in  which  the  int.ens'-  field  is  detuned  from 
the  center  of  its  broad  Doppler  profile.  Considerable  simplification  then 
results,  and  the  standing-wa\  response  may  be  analyzed  in  terms  of 
pairs  of  travelling- wave  fields  interacting  with  the  respective  transitions 
of  the  three  level  system.  As  will  become  evident,  the  relative  propaga¬ 
tion  direction  of  p:  one  and  saturating  field  components  is  of  crucial 
importance:  fields  propaguiing  in  the  same  direction  lead  to  a  probe- field 
line  shape  which  is  strikingly  different  from  that  due  to  oppositely-pro¬ 
pagating  fields. 

The  atoms  of  a  Doppler-broadened  gas  may  traverse  many  wave¬ 
lengths  of  the  applied  fields  before  decaying.  We  adopt  here  the  simple 
picture  of  atomic  motion  in  which  an  atom  produced  in  a  particular  state 
travels  undeflected  with  constant  velocity  as  it  decays.  The  calculation 
of  induced  emission  may  be  divided  into  two  stages:  The  response  is  first 
obtained  for  a  band  of  three  level  atoms  within  a  narrow  range  of  axial 
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velocities  interacting  with  the  pp|i<d  travel  ling- wave  fields;  this  f|uanlity 
is  then  summed  over  the  entire  distribution  of  velocities,  thus  obtaining 
the  complete  emission  profile.  The  first  stage,  the  calculation  of  the 
ensemble-averaged  travelling- wave  response,  may  be  carried  out  in 
several  ways.  The  induced  dipole  moment  of  an  atom  produced  in  a 
given  state  may  be  calculated  from  the  Svhr’o’dinger  equation.  Equiv¬ 
alently,  one  may  calculate  the  rate  at  which  an  atom  produced  in  a  given 
state  makes  transitions  to  other  states.  In  either  ease  the  response  of 
the  entire  velocity  ensemble  is  obtained  by  averaging  the  quantity  cal¬ 
culated  over  all  initial  conditions.  The  transition  rate  approach,  which 
has  the  important  advantage  of  identifying  the  various  radiative  processes 
by  which  an  atom  emits  and  absorbs  photons,  is  presented  in  Section  1IA. 
The  connection  between  transition  rates  and  induced  dipole  moments  is 
examined  in  Appendix  B.  An  alternate  derivation  of  the  line  shape  using 
the  ensemble-averaged  density  matrix  equations  of  motion,  and  related 
discussions,  are  presented  in  Appendix  C. 

A.  Ensemble- Averaged  Response:  Transition  Rate  Approach 

The  resonant  interaction  of  two  monochromatic  fields,  P’^  and 
with  a  three  level  system  was  treated  in  Ref.  10  for  eases  where  Doppler 
broadening  is  negligible  (e.  g. ,  the  microwave  region)  and  the  decay  rales, 
considered  equal  for  dll  three  levels,  were  assumed  to  result  from  hard 
collisions.  The  perturbation  method  consisted  of  first  obtaining  a  closed- 
form  solution  to  the  Schrodinger  equation  for  E^-  0  and  E2  arbitrary,  and 
then  using  this  result  to  generate  a  solution  valid  to  first  order  in  E ^ .  Tht 


I 
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present  section  is,  in  part,  a  generalization  of  that  method  valid  for  levels 
with  differing  decay  rates. 

Ti  e  three  level  systems  to  he  studied  are  of  the  type  shown  in 
Fig.  2(a),  (b),  (c).  Level  0,  the  common  level,  is  coupled  to  levels  1  and 
2  bv  e  lectric  dipole  matrix  elements  and  r'('sPt‘fdively.  Denote 

the  energy  of  level  j  by  hW.,  and  let  I  W . -  Wq|  -  w.;  and  fall  in  the 
optical- infrared  region  we  shall  assume  (u^  -  Wg)  is  large  compared  to 
the  natural  linewidths.  In  the  present  Section  we  consider  the  "inverted- 
V”  level  configuration,  in  which  level  0  lies  highest  [Fig.  2(a)].  The 
treatment  is  easily  extended  to  cases  in  which  level  0  lies  below  either  or 
both  of  levels  1  and  2.  This  extension  is  discussed  in  Appendix  C. 


The  system  interacts  with  a  strongly  saturating  field  F^(z,t:«)  at 
a  frequency  close  to  w,^.  Tie  resonance  at  is  probed  by  the  weak 
field  Kj  (x,  t)  at  variable  frequency  To  allow  for  both  possible  relative 

propagation  directions,  Fj(/,t)  is  taken  co  be  travelling  in  the  positive 
{+•/.)  direction,  while  K.^(y.,\  \  ()  may  propagate  in  either  positive  (c  -  L  1)  or 
negative  (c  -  1)  dirertions.  Speeifieally, 


and 
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Consider  an  ensemble  of  atoms  at  given  position  and  time,  and 
moving  with  axial  velocities  in  the  narrow  interval  between  v  and  v+dv. 

In  a  coordinate  system  in  which  the  ensemble  is  at  rest  the  incident  fields 
appear  as  E '  (z ,  t )  and  El,(z,t;€),  with  E1.  identical  to  Eqs.  (1)  except  that 
due  to  the  l)oppl<  r  effect, 

ilj-  =  Sij  -  k,  v  , 

n2  =  n2‘  ek2V  ;  (2) 

k.~  k*.  =  n'./c 
J  J  J 

Note  that  in  this  section  z  and  t  always  refer  to  the  coordinates  in  the 
moving  frame.  The  total  Hamiltonian  for  the  system  is 


H  =  H0  +  V  (t) 


(3) 


where  HQ  is  the  Hamiltonian  of  an  isolated  atom  with  stationary  states 
of  energy  fi\V d: 


HQ^j  =  MV.tf/.  ,  j-- 0,1,2 


(4) 


The  ensemble  is  coupled  to  the  applied  fields  by  the  interaction  Hamiltonian 


V  (t)  -  -  mE'  (z.t) 


(5) 


in  which  p  is  the  electric  dipole  operator  and  E*  -  Ej  +  E^.  The  time 
evolution  of  a  particular  member  of  the  ensemble  is  determined  by  its 
wave  function 't,  which  may  be  expanded  in  terms  of  the  stationary  states 


V 
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*  (R,t)  =  N  e  (t)  0  (R)  ,  (6) 

L-t  J  .> 

j  =  o 

v.  ith  tj(t)  the  probability  amplitude  of  level  j.  The  equations  of  motion 
for  the  e.’s  may  be  obtained  in  the  usual  manner  by  inserting  Eqs.  (3)- 
(6)  into  the  time-dependent  Schrodinger  equation  U#  -  i  fid'P/St, 
multiplying  on  the  left  by  ,  and  integrating  over  all  space  (R).  Going 
over  to  the  interaction  picutre  by  means  of  the  transformation 

iW  t 

Cj(t)  =  c.(t)e  J  (7) 

one  obtains: 

6.  =  )  a. .  c.  on 

i  L  iJ  3 

.1 


in  vdiieh 


a. 

i.l 


-  V . .  e 

h  ij 


i(\Yj  -  Wj)t 


(9) 


vir'“ijE'  (F’t> 


(10) 


(In  the  above  equations  R  indicates  the  electron  coordinates  in  the  atom’s 
reference  frame  and  r  is  the  position  vector  of  the*  atom's  center  of  mass 
in  the  moving  frame.  ) 

Equation  (8)  describes  the  evolution  of  an  undamped  atom.  The 
effects  of  radiative  decay  may  be  included21  *by  modifying  (8): 


=  /  (a.  -  —  y  6.  )  c . 
1  L  'J  2  J  t.l  J 


(ID 


Note  that  in  the  absence  of  applied  fields  (11)  leads  to  exponential  decay  of 


level  j  with  decay  rate  yr. 


|c.(t-t»)l2  =  |C  (t')|2  €•  J 
3  J 


excitation  arising  from  incoherent  processes.  The  latter  are  responsible 

for  populating  the  levels  at  rates  which  may  be  assumed  to  be  independent 

(221 

of  the  applied  fields,  accordingly,  background  atoms  are  produced  in  level 
k  at  a  rate  n^y^,  where  n^  is  the  number  of  atoms  with  velocity  component  v 
in  level  k  in  the  absence  of  applied  fields. 
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Wo  adopt  the  notation  Cj(t;t0,  k)  to  indicate  a  solution  to  (11)  for  an  atom 
produced  in  level  k  at  The  appropriate  boundary  conditions  are: 


c  .(t  :t  ,  k)  -  6.  .e^j  ’ 
j  o  o’  kj  J 


(16) 


in  which  the  (p.'s  are  unimportant  arbitrary  phase  factors.  The  transition  probability 
3 

2 

to  another  level,  level  j,  at  subsequent  time  t  is  i  o t : tf>,  k) !  The  prob¬ 
ability  of  an  atom  in  level  j  decaying  ir:  the  interval  between  t  and  t  +  dt  is 
y  .dt.  Thus,  for  an  atom  produced  in  level  k  at  t  ,  the  transition  rate  to  level 

/is:(24) 


(17) 


Considering  that  the  rate  of  production  of  atoms  in  level  k  is  the 

ensemble-averaged  transition  rate  between  levels  k  and  j  is: 
t( 

nk  7kTj  \  kj  (t;  t(),  k)l2  dtf)  nk.Ikj  ,  (18) 

-  S 

which  defines  .1,  the  ensemble-averaged  k— *  j  transition  rate  per  atom 
KJ 

produced  in  level  k. 

In  order  to  estimate  the  net  emitted  power  induced  by  E^,  all 
events  in  which  a  photon  is  emitted  or  absorbed  at  n'j  must  be  considered; 
An  atom  which  is  produced  in  level  2  and  subsequently  decays  from  level 
1  at  rate  must,  by  necessity  of  energy  conservation,  exchange  two 
photons  with  the  applied  fields,  a  photon  absorbed  at  an(^  a  photon 
emitted  at  [  Fig.  2  (a)]  .  In  contrast,  an  atom  which  is  produced  in 
level  0  and  subsequently  decays  from  level  1  at  rate  Jqj  must  emit  a 
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single  photon  at  12^.  An  atom  produced  in  level  i  absorbs  photons  at 
li'j  by  the  reverse  processes,  namely  single-quantum  transitions  from 
levels  1  to  0  and  double- quantum  transitions  from  levels  1  to  2.  The  net 
rate  of  emission  of  photons  at  is  therefore 

($  =  n2J21  +  n0J01  '  nl  ^10 +  J12^  ’  (19a) 

and  the  corresponding  power  emitted  is: 

iv(«i)  =  •  (19b) 

Note  that  in  the  limit  of  K' —  0,  .J,^  and  should  vanish,  and  J^j  and 
should  reduce  to  the  usual  transition  rates  for  a  two  level  system. 

In  the  formulation  presented  here  the  elementary  act  of  transition 

from  level  k  at  t0  to  1<  vel  j  at  t  is  described  in  terms  of  its  corresponding 

2 

transition  probability,  |c.(t:t0,  k) I  .  It  is  evident  on  ”ery  general  grounds 
that  the  reverse  act  of  transition  is  equally  probable; 


|e.  (t;l0,  k)r  =  !ck(t;t0,  j)r  .  (20) 

This  point  is  discussed  further  in  Appendix  A.  It  immediately  follows  that 
Jki  ~  and  (10a)  simplifies  to 

(n2-  nj)  J21  +  (nQ-  nj)  JQ1  ,  (*21) 


a  convenient  form  for  calculation. 

It  is  v.orth  noting  that  F.q.  (19)  utilizes  the  fact  that  the  atomic 
ensemble  on  the  average  exchanges  electromagnetic  energy  with  Kj  in 
units  of  This  is  interesting,  since  throughout  our  formulation  we 

have  dealt,  with  classical  fields  exclusively.  In  fact,  in  calculating  the 
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emitted  power  it  is  not  necessary  to  explicitly  introduce  the  concept  of 

(3) 

field  quanta;  instead,  one  may  evaluate  the*  expectation  value  of  the 

('  j|:  — *  3 

induced  dipole  moment,  \  'k  e  R  d‘  R.  Then  the  ensemble-averaged 
polarization  at  ft.  for  atoms  produced  in  level  k  is 

t 

Pjk(t)  =  nk  \  \  2Re  [Mj0e::(t;tOJk)e0(t;to,k)jdto  ,  (22) 

-  yi 

and  the  net  power  emitted  is 


Iv(S!j)=‘<b^k<t,Ej  <Z,t)  ’tin- average 


(23) 


In  Appendix  B,  using  density  matrix  notation,  it  is  shown  that  Eq.  (23)  leads 
to  I  (fl'j)  =  hft!^  ,  Eq.  (1 9b),  where ^  is  identical  to  Eq.  (19a).  This  is 
accomplished  by  calculating  the  power  emitted  or  absorbed  by  atoms  pro¬ 
duced  in  particular  levels,  starting  from  Eq.  (22).  For  instance,  in  this 
way  it  is  shown  that  for  atoms  produced  in  level  2  [Fig.  2(a)],  the  emitted 
power  at  ft'^  is  given  by  hft’ n^.J^j,  while  for  atoms  produced  in  level  0  it 
is  given  by  bft'j  HqJqj. 

We  proceed  to  calculate  the  emitted  power  in  terms  of  transition 
rates,  Eq.  (21).  The  explicit  form  of  the  equations  of  motion  (14)  are: 


.  i  6. t  i 6  q  t 

dQ  =  a*  e  dj  +  i)3  e  d 2 

-i6 j t 

dj  =  iae  dQ  , 


*  "i62t 
=  i|3*e 


(24a) 

(24b) 


(24c) 
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in  which 


a  =  Hj  0Aj/t> 


(25) 


(3-^02  A2 /fl 


and  i6^  =  i(w^-  ^ )  +  —  (Tq-  7^).  In  writing  (24)  anti-rcsonant  terms  have 

2 

been  neglected,  since  their  influence  is  negligible  when  <pE>  <<  t> “ 
From  Eqs.  (13)  and  (16),  the  appropriate  boundary  conditions  are  seen  to 
be 


(,j  ttoi'o*  k)  6kj 


-  y.  t,  +  i <p. 
2  .1  °  *.l 


(26) 


A  solution  to  this  set  of  equations  for  all  values  of  E9  and  for  weak 
Ej  may  be  obtained  by  means  of  a  simple  perturbation  technique. This 
procedure,  however,  involves  lengthy  expressions.  To  present  the  im¬ 
portant  steps  in  a  concise  manner,  a  number  of  symbols  are  introduced. 
For  the  convenience  of  the  reader,  these  are  collected  in  Table  I.  Addi- 
onal  details  regarding  the  perturbation  technique  will  be  found  in 
Appendix  A. 

Consider  an  atom  produced  in  either  level  0  or  2:  In  the  absence 
of  coupling  through  a  (i.  e.  a  -  0),  (24b)  yields  dj  =  0,  while  (24a)  and  (24c) 
reduce  to  the  equations  of  motion  of  a  damped,  two  level  system.  The 
solution  is  straightforward:  for  an  atom  initial!;,  in  level  2, 


i'W  i 
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<(Nt:t  ,  2)  = 


i  "o  t 


V<‘ 


iq_T 

)  , 

(27a) 

-  i  d+  T 

)  . 

(27b) 

,u  /  .  lirj  (t, tn,  2)  ,  -iO_T 

( 2  t;  tQ‘  2)  r  2-“-°  (  %  -q_(1 

In  these  equations  the  superscript  u  designates  the  parameters  of  the  un- 

„  /  9  2 1  l 

coupled  system  and  T  =  t  -  tQ,  s  =  V°^  +  4  1/3  i  and  q+  ~  (  o9  +  s).  By  virtue 

of  symmetry,  d)?(t;t  ,  0)  may  be  obtained  from  d!?(t;t  ,  2)  bv  replacing 
d2(to;to,  2)  by  d0(xo;to,0)  and  q+  by  -q_.  For  a  i  0,  d}  no  longer  vanishes: 
an  approximate  expression,  complete  to  lowest  order  in  a,  may  be  obtained 
by  integrating  (24b)  with  dp  replaced  by  d^: 


C  'i6^v  „ 

dj(t;to,k)  a  ia  \  e  djj  <tf  :tQ,  k)  dt1 ;  k  =  0,  2  . 


(28) 


For  tn  atom  initially  in  level  2,  substitution  of  (27a)  into  (28),  and  using 
(13),  yields; 


C1  lo 


;  t  ,2)=  9Lii  c"2  ?1  Te  {(^21’  ",21)to"V>2  }je 


iq^l  in  T 

_ zl.j ; _ ii\ 


9  + 


(29) 


in  which  Sl9,  =  fig'  u2l~  tu2~  u\’  an(*  9+ "  Oh  _  ^ i  •  Similarly,  (t;  tf),  0) 
may  be  obtained  by  substituting  dj.1  (t;tQ,  0)  into  (28);one  finds: 

t  >0)=g.e"t  q+  [ e  1  r|-T- lj  q_  [cir'+T-l)l  , 


c  i  (t  ■  *o’  ^  s 


n_ 


Hh 


,  (30) 


r  Symbols  Used  ?n  Calculation  of  Transition  Rates,  Eqs.  (31)  and  (32) 


°21  =  S22  ‘  n'l 

a  =  M10 Ai^ 

W21  r  w2' 

ft  =  HQ2  a2  /h 

i  7-  t 

0.  k)  -  r  .  (t;  tQ,  k)  0  1  J 

T  =  l-t0 

i(wk"  “k)  +  |  (yQ~  V 

q  +  "  ^  (62  ^  s) 

2+4  ;3i2 

ri+  = 

-  22  - 


ith  (18),  J  and  J  may  be  obtained  by  integrating  I  tH  , 
01  21 


In  accordance  with (18),  J^ana  .  - 

obtained  from  Eqs.  (29)  and  (30),  over  t0,  one 

M2!(3!2  WfJ  _J_  |  i  +_ - - - — - — - ! 

j21  =  7i72  jTfT-  i  |n+!2  1  7i  Vi(n+_T1+  7i‘l71  + 


i  i  + 


+  |  .2  17!  7i-i(n-9i!>  'Yl'19+- 


1+ _ J 


77  Vi-HrU-TU*)  7l'iTl- 


1  i  \  . 

71+iq+*J 


2  2 

•T01 '  77?  lie  'hJ7  '  7  ‘ n+ 


+  17^1  iFT^rXn.^)  n-in-J 


A-  I 

1-  IT!_  J 


The  const, tuen.  terms  of  (3U  and  of  (32)  may  oe  con,,, 
a  straightforward  (though  lengthy)'  manner  to  obtain: 


be  combined  ~nd  simplified  in 


720 

0  9  i  L2  ” 2  77T  \  . 

J21=  2  531  Im^  Ib  ’ 

L2-2^R  ^ 

2  ; R  ,o,2  2  ^0  l 

Jor-21"1  Im  -B  p  "AB 
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Hero, 


A  "  I  L0  1 2  +  — .g?  \'i  ;2 


y0y2 


n  =  ~K  Lj  +  j(i  j 


Li  '  Ai+  ‘I'lo 


h  = ~ A9+17 


2  r20 


H= 


21 


A'.  =  Si'  -  u 
J  J  J 


and 


ij 


i(vv 


As  exported,  for  |0|2~  0,  J2J  approaches  zero  while  J„.  reduces  to 

2,2  U1 

2  *io'  *  *Jl  ^  usua*  expression  for  the  0—  1  single-quantum 

transition  rate.  Accordingly,  ‘.he  laser  field  nanifests  itself  through 
two  separate  radiative  acts:  1),  it  gives  rise  to  a  double -quantum  tran¬ 
sition  rate  J2l;  2),  it  modifies  the  single-quantum  transition  rate  JQ1. 
As  discussed  in  Section  IIIA,  below,  for  tuned  close  to  u  and  n' 
cluSe  to  the  magnitude  of  J,,,  is  of  the  same  order  as  the  H- 
cleponden,  portion  of  J#r  Atd,.tuned  frequencies,  however,  J2[  mav 
be  considerably  larger  than  J*.  fact.  „  rk.,uned  frequeneie 
reduces  to  the  familiar  expression ,z5>  for  Raman  transitions  be- 
tw  een  levels  2  and  1.  [See  Section  IIIA  for  further  details.  ] 


(35a) 

(35b) 
(36  p.) 
(36b) 

(36  c) 

(37) 

(38) 


ics  J 
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Having  obtained  explicit  expressions  for  Jgj  and  the  net  rate 
of  emission  at  £2^  in  the  moving  frame  follows  from  Eq.  (21): 


^  =  2  I  a  I  2  Im£)  (v,  €  ) 

^20 

JD(v,  £  )  =  K-n^t  (n2  -  nQ)  1/3  1 2  --  - 


As  in  previous  notation,  the  variable  t  specifies  whether  Ej  and  Eg  propagate 
in  the  same  direction  (c  =  +  1)  or  in  opposite  directions  (c  -  -  1).  In  the 
laboratory  rest  frame  photons  are  emitted  with  energy  Thus,  in  that 

frame  the  contribution  to  the  emitted  power  from  atoms  moving  with  axial 
velocity  v  is 

Ti  Q  2  |a  1 2  Im  J)(v,£  )  .  <41> 

Setting  all  the  equal,  (41)  reduces  identically  toexpressions  obtained  in  Ref.10, 

where  hard  collisions  were  introduced  as  the  mechanism  of  decay.  The 

(26) 

detailed  features  of  that  line  shape  have  been  fully  verified  m  the  micro- 
wave  region  where  velocity  broadening  is  negligible.  The  reader  is  referred 
to  Refs.  10  and  26  for  additional  details. 

In  the  formulation  presented  above  the  saturation  effects  manifest 
themselves  as  nonlinear  intensity  dependences  of  the  transition  rales.  The 
level  populations,  n . ,  enter  into  the  expression  for  the  net  emission  rate 
through  the  background  excitation  rates  n.y.,  which  are  assumed  independ¬ 
ent  of  the  applied  fields. Nevertheless,  the  ensemble  averages  of  these 


populations,  tv,  do  depend  on  the  intensities  of  the  applied  fields.  These 
average  populations  may  be  computed  in  a  straightforward  manner  by  con¬ 
sidering  the  average  number  of  atoms  found  in  level  j  due  to  production  of 
atoms  in  all  possible  initial  levels  k: 
t 

■a  i  \ 

(42) 


"j  '  }  "kl'k  j  |cj  <*;  tQ,  k>l2  0t„  1  N  „kJk.  . 


J  k 


The  right  hand  expression  immediately  follows  from  the  definition  of  J^, 
Eq.  (18).  It  should  be  noted,  however,  that  the  average  level  populations 
n.  do  not  enter  explicitly  into  the  expressions  for  the  power  emitted  or 
absorbed  at  and  That  these  quantities  are  related  to  the  net  rate 

of  photon  emission  may  be  seen  as  follows:  It  is  shown  in  Appendix  B 
[Eq.  (B  1 9b )]  that  the  transition  rates  are  connected  by  the  relation. 


Yi  '  ^  'V 

k 


(43) 


This  equation  is  merely  a  restatement  of  the  ,-teady-. state  condition:  Rate 
of  production  into  level  j  equals  net  rate  of  decay  from  all  levels.  Equa¬ 
tions  (42)  and  (43)  may  be  combined  to  yield 

^  (nkJkj-  njV =  ®r  •  (44) 

k 

For  j  -  1,  the  left  hand  side  is  Eq.  (19a).  Accordingly,  this  ex¬ 

pression  may  be  interpreted  as  the  steady-state  form  of  a  rate  equation 
for  the  population  of  level  1 , 
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nt  =  0  =  (lij  -  n1)71  +  ^  ,  (45) 

underscoring  the  fact  that  emission  and  absorption  of  photons  must  be 
accompanied  by  atomic  decay. 

As  pointed  out  earlier,  an  equivalent  way  of  calculating  the  emitted 
power  at  ST2 is  to  solve  the  ensemble-averaged  density  matrix  equations  of 
motion  in  the  steady  state  and  obtain  the  induced  polarization.  In  tht  latter 
approach  one  does  not  generally  distinguish  between  single-quantum  and 
double- quantum  processes.  Also,  the  solution  of  the  density  matrix  equations 
involves  algebraic  manipulations  quite  different  from  those  used  above. 

This  approach  is  presented  in  Appendix  C.  Tv'e  method  of  solution  is  more 
or  less  standard,  except  that  in  the  perturbation  approach  employed  one  of 
the  applied  fields  may  be  taken  as  arbitrarily  large.  The  final  results 
are,  of  course,  in  complete  agreement  with  Eq.  (39). 

B.  Doppler-Broadened  Response 

Having  obtained  the  power  •  mitted  bv  the  atoms  moving  with  axial 
velocity  v,  it  is  now  necessary  to  sum  over  the  entire  distribution  of  ax:  >l 
velocities,  thus  obtaining  I  (f.2 ^ ;  c ),  the  total  emitted  power  at 

y. 

I  («i;e  )  =  2h«1  !o|2  Im  (v,  e)dv  .  (4c  / 

In  carrying  out  the  velocity  averages  it  is  convenient  to  introduce  G.(v), 
the  velocity  distribution  of  atoms  in  level  j: 

G j  (v)  d  \  =  1 


n .  -  N  G  (v) 
.1  J  J 


(47) 
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with  Nj  the  total  number  of  atoms  in  level  j.  To  an  incident  light  beam  of 

propagation  constant  k  the  frequency  breadth  associated  with  may  be 

characterized  by  a  "Doppler  width"  ku.,  with  u.  as  the  most  probable  speed 

of  G..  Where  the  velocity  distributions  are  thermalized,  the  G.(v)'s  are 
1  '  J 

Maxwellian.  In  some  applications  of  Laser-  Induced  Doppler  Line  Nar¬ 
rowing,  however,  the  •  elocity  distributions  are  non-thermal  and  may  even 
deviate  considerably  from  Gaussian  form.  (See,  in  particular,  discussions 
of  the  atomic  oxvgen  laser,  Hofs.  3,  7,  and  20.  )  In  general,  the  integrals 
involved  in  obtaining  (46)  depend  upon  the  specific  form  of  the  Gds. 
Nevertheless,  because  in  the  present  case  we  are  considering  the  fully 
Doppler-broadened  limit  y/ku  «  1  [u  and  y  characterized  the  magnitudes 

of  the  most  probable  speeds  (u  ’s)  and  the  natural  widths  (y.,s),  respectively) 

J  J 

the  resonant  behavior  of  I(S2j;e)  becomes  largely  independent  of  the  Gds 
and  it  is  possible  to  perform  the  averages  without  recourse  to  their  specific 
form. 

The  velocity  averages  are  carried  out  in  detail  in  Appendix  D. 

Briefly,  the  velocity  dependence  ofjD(v:c),  which  enters  through  Lj (O'j) and 
Lgffig).  Kfl  (36a,b), and  Rfnj,  Eq.  (36c),  as  well  as  n-,  Eq.  (47),  is  rather 
complicated  in  its  present  form;  therefore,  as  a  first  step, $)  is  rewritten 
in  a  form  consisting  of  verms  having  velocity  dependence  of  the  type 


H7  (v)  = 


G.  (v) 


(u  +  iy)  +  k  v 


(w •  +  iy'  )  +  k'v 


V 


(48) 


in  which  w,  u>',y  and  y'  are  real  and  y  and  y’>  0.  As  shown  in  Appendix  D, 
in  the  Doppler  limit  y/  ku  «  1, 


A  .  “  S2 .  -  w . 

J  .1  .1 


(52) 


Q  is  the  factor  by  which  the  saturated  levels  of  the  0-2  transition  art; 
broadened: 


and  t  specifies  whether  and  E9  propagate  in  the  same  direction  (c  +1) 
or  in  opposite  directions  («  =  -  1 ). 


Equation  t50)  predicts  a  sharp  Lorentzian  decrease  (for  N.-,) 

superimposed  upon  the  broad  gain  profile,  an  effect,  due  to  the  nonlinear 

coupling  of  Ej(z.t)  and  Eg(z,t;e).  The  center  frequency  and  width  of  this 

resonant  decrease  are  dependent  upon  the  relative  propagation  directions 

of  the  applied  fields.  When  E.  and  Kg  propagate  in  the  same  direction 

(«  =  +  1  case)  the  change  signal  is  detuned  front  Uj,  the  center  of  the  0-1 

<a. 

Doppler  profile,  bv  an  amount  A.  =  — -1-  A0  and  is  of  width 

1  u2  1 

rN',dl  ■  <54> 

On  the  other  hand,  when  Ej  and  Eg  propagate  in  opposite  directions 

(«  =  -  1  _ase)  the  change  signal  appears  at  a  frequency  detuned  from  by 

u. 

A,  -  — *-  A„  and  is  of  width 
1  2 

V  Ti  + 1  jc>  •  <55> 

which  is  broader  than  r^.  by  2  7qQ,  twice  the  saturated  width  of  the  middle 
level.  The  frequencies  at  which  the  c  -  +  1  and  t  =  -  1  change  signals  occur 
have  the  following  significance.  Due  to  the  Doppler  effect,  each  applied 
field  couples  resonantly  only  to  atoms  within  a  narrow  band  of  velocities. 

P  ferrirg  to  Eq.  (2),  it  can  be  seen  that  E ^  resonates  with  atoms  of 
velocity  mar  ,  given  by  vj  ~u\-  Similarly,  Eg  resonates  with 

atoms  of  velocity  near  v(>,  given  by  Ilg-  ek^Vg-  Ug.  [n  general,  the 
velocity  bands  centered  at  and  v,;  are  distinct  and  do  not  overlap.  How¬ 
ever,  at  particular  frequencies  of  ih"  applied  fields,  the  two  bands  will 
merge  into  a  single  one.  Equating  \j  and  Vg  for  t  -  +  1,  this  condition  is 
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srf  n  to  be 

Ai  =  €rrA2  *  <f)6> 

which  is  the  location  of  the  center  frequency  of  the  corresponding  change 
signal. 

Evidently,  the  applied  fields  can  couple  to  each  other  most  effectively 
when  they  resonate  with  the  respective  transitions  of  atoms  within  a  narrow 
band  of  velocities.  Complete  substantiation  of  this  point  follows  from  a  de¬ 
tailed  inspection  of  the  frequency  behavior  of  the  rest  frame  response- - 
determined  by  the  corresponding  transition  rates--  tor  a  particular  atomic 
velocity  band.  A  full  discussion  is  deferred  to  Section  III  A .  Several  re¬ 
marks,  however,  may  be  useful  at  this  time.  A  glance  at  Eqs.  (33)  and  (34) 
reveals  that  Jgj  and  are  rather  complicated  functions  of  and  d!>. 
Furthermore,  the  frequency  characteristics  of  JQ1  are  quite  different  from 
those  of  J9j.  In  the  fully  Doppler-broadened  limit,  however,  a  number  of 
simplifying  cancellations  occur  and  many  details  of  the  responses  of  in¬ 
dividual  velocity  bands  average  out.  In  fact,  in  this  limit  the  velocity- 
averaged  2—  1  transition  rate  becomes  equal  to  the  negative  of  the  ft- 
dependent  portion  of  the  0--  1  transition  rate,  a  fact  responsible  for  the 
particularly  simple  form  of  the  final  expression  for  the  emitted  power, 

Eq.  (50). 

Despite  the  fact  that  the  €  r  +  1  and  c  =  -  1  change  signals  are  symme¬ 


trically  located  about  the  0  -  1  center  frequency,  the  corresponding  widths 
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differ.  This  asymmetry  results  from  the  resonant  behavior  which  enters 
,J21  and  through  the  quantity  R  =  [(fij-  ft!,)  -  1  ~  iY2l  = 

l  12”  '  YS 1 1  "  v  ^i"  *  k2).  Note  that  in  contrast  to  Lj  and  I^.  the 

velocity  dependence  of  R  differs  strongly  for  c  =  +  1  and  «  =  -  1,  particularly 
if  Wj  and  are  comparable.  For  example,  in  the  important  case  in  which 
I  <^2~  f  v/c  «  Y21  ’  R  {«  =  +  1  }  becomes  essentially  velocity  independent, 
while  R  {«  =  -  1}  remains  strongly  velocity  dependent.  Bearing  this  fact  in 
mind,  the  different  widths  of  1'^  and  Tp  can  be  understood  by  inspecting 
(33)  and  (34).  See  Section  HI  A  for  further  details. 

As  discussed  in  Appendices  C  and  D,  the  preceeding  remarks  also 
apply  to  the  other  level  configurations  of  interest.  For  example,  the  ex¬ 
pression  for  emitted  power  for  the  "V"  configuration  [Fig.  '2(b)]  ,  which 
is  the  inverted  version  of  the  level  scheme  treated  in  this  section,  is 
given  by  the  negative  of  Eq.  (50),  as  one  might  expect.  Equation  (50)  also 
describes  the  power  emitted  by  the  cascade  configuration  [Fig.  2(c)]  when 
the  factor  £  Q  appearing  in  the  imaginary  part  of  the  denominator  is  changed 
to  -  €  Q.  This  has  the  interesting  consequence  of  interchanging  the  positions 
of  the  broad  and  narrow  change  signals  (for  a  given  value  of  fl9),  a 
"geometrical"  effect  arising  from  the  different  way  in  which  the  fields 
co1. pie  to  a  cascade  system. 

C.  Extension  to  Standing- Wave  Fields 


The  extension  of  the  above  considerations  to  cases  where  one  or 
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both  of  the  applied  fields  are  in  the  form  of  standing  waves  is  straight¬ 
forward.  As  discussed  above,  nonlinear  coupling  occurs  when  two  (or  more) 
travelling- wave  fields,  at  least  one  of  which  can  saturate  its  transition,  are 
Doppler-shifted  into  resonance  with  the  >ame  narrow  velocity  band  of  atoms. 
As  noted  in  the  introduction,  each  standing- wave  field  may  be  decomposed  into 
travelling-wave  components  of  equal  amplitude  propagating  in  opposite 
directions.  As  long  a  the  strong  field  is  detuned  from  the  center  of  its 
broad  Doppler  profile  (i.  e.,  '  -  oj  !  (;,)mP°nents  resonate  with 

atoms  in  distinct  bands  symmetrically  located  about  the  center  of  the  velocity 


distribution.  (In  this  respect  the  tuning  of  a  weak  standing-wave  field  is 
unimportant,  since  its  travelling- wave  components  cannot  couple  to  one  an¬ 
other.  )  Then  the  nonlinear  coupling  can  occur  only  when  the  frequencies  of 
the  applied  fields  are  such  that  one  of  the  travelling-wave  components  of  the 
intense  field  and  one  of  the  weak  field  are  Doppler- shifted  into  resonance  with 
the  same  hand  of  utems.  One  such  possibility  couples  together  travcdling 
waves  propagating  in  the;  same  directions  arid  will  be-  recognized  as  the 
c  *  1  condition  of  Kq.  (56).  In  aodilion,  there  is  the  possibility  of 
coupling  between  oppositely-propagating  travelling  waves,  which  is  just 
the  *  -  -  i  condition  of  Eq.  (56).  Consequently,  the  broad  and  narrow 
resonances  appear  simultaneous  y,  symmetrically  located  about  the 
0-  l  center  frequency. 

To  illustrate  the  preceding  remarks,  consider  an  intense  standing- 
wave  field  of  amplitude  2  and  a  weak  travelling-wave  field  of  amplitude 
resonating  with  the  re  spec  tive  transitions  of  a  Doppler-broadened  three 
lev*,  i  system.  For  the  present,  assume  to  be  detuned  from  u,,.  This 


situation  is;  directly  applicable  to  the  example  presented  in  the  introduction 


where  probes  the  linoshape  of  an  optical  transition  sharing  a 

common  level  with  an  oscillating  gas  laser  transition  (Fig.  1). 

I  or  concreteness,  suppose  the  velocity  distributions  Gj(v)[Hq.  (47)  |  are  all 
Maxwellian  at  temperature  T:  In  the  absence  of  the  laser  field  (Kq  =  0), 
^(i‘2j),  the  spectrum  of  power  emitted  at  the*  probe  frequency,  is  just: 

A  (N0-N,)g&>lf  n~(|Vi)  .  (57) 

v  1  1  u  kjU/tyjf 


which  defines  ■Jfy  (ft, ),  the  usuai  expression  for  power  emitted  by  the  Doppler- 

broadened  0-1  transition  induced  by  weak  travelling-wave  field  E. ;  here, 

1  '2 

u  -  (2k  T  M)  is  the  most  probable  speed,  M  is  the  atomic  mass,  and  k 
is  Boltzmann's  constant.  In  the  presence  of  the  laser  field,  ^(ft^)  is 
strikingly  modifif  d:  from  the  c  -  +1  and  c  =  -  1  eases  of  Eq.  (50)  one  obtains. 


^  (nt)  $  ('-2, )'  1  +  \ 


VN1 


VqQ 


Im  J 


1 


(A  A,)+i^B  (A  -a 

1  Wo  J  O  1 


W2 


V 


+  i 


■h 


so  1  hat  narrow  and  broad  change  signals  identical  to  those  described  above 
for  the  travelling-wave  ease  appear  simultaneously  on  opposite  sides  ofwj. 

In  a  number  of  important  applications  T1rr7’2  <<;  Tg-  causin8  ^  and  %  to 
differ  enormously. 

Equation  (a8),  which  was  originally  presented  in  Ref.  7  ,  has  been 
written  in  a  form  valid  for  both  case. -me  (n  ~  -  1)  and  "inverted-V"  (a  =  +  1) 
coupling  configurations,  Figs.  2(c)  and  2(a),  respectively.  The  corresponding 
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expression  for  the  "V"  level  srheme,  Fig.  2(b),  is  given  by  the  negative 
of  (58)  with  7  =  +  l.  These  statements  are  based  on  the  extension  of  Eq.  (50) 
to  the  other  coupling  schemes,  as  outlined  in  Appendix  D.  See  also  the 
remarks  at  the  close  of  Section  IIB. 

In  the  case  of  weakly  saturating  laser  field,  reduces  to 

CO 

2'y10  +  wi  720‘  V  and  PB  to  2^710+[TL  y2<?'  which  a8rees  with  exnres- 

^  ( r  6 ) 

sions  previously  obtained  by  other  methods.  '  It  is  easily  seen  that 
in  that  limit,  by  symmetry,  the  case  of  u>2>  o.  may  be  obtained  by  simply 
interchanging  and  in  (58). 

The  interpretation  of  (58)  as  it  applies  to  the  spontaneous  emission 
version  is  presented  in  Section  IFB. 

Equation  (58)  is  valid  as  long  us  the  intense  standing-wave  field  is 
detuned  from  its  atomic  center  frequency.  When  both  fields  are  tuned  to 
the  centers  of  their  respective  gain  profiles  all  of  the  travelling-wave  com¬ 
ponents  can  couple  to  the  same  atomic  velocity  band  (namely,  the  one  with 
negligible  component  of  velocity  along  ihe  z-axis).  Consequently,  the  above 

considerations,  which  are  based  on  the  coupling  of  pairs  of  travelling  waves, 

(  2  ) 

becomes  inadequate.  A  standing-wave  analysis  which  specifically  includes 

this  possibility  has  described  the  third -order  interaction  of  E  and  E  for  the 

1  2 

special  'ase  in  which  levels  1  and  2  of  the  "V"  configuration  (Fig.  2(b)) 
are  assumed  closely  spaced  (I  j  v/c  «  y).  [Note  that  the  analysis 

of  the  cascade  configuration.  Fig.  2(c),  in  which  the  middle  level  lies 
about  half-way  between  the  other  two  levels,  follows  identically.  )  The 
analysis'  ,  based  on  a  density  matrix  calculation,  shows  th.-it  as  is 


lllUMi 
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tuned  close  to  the  center  irequency  of  the-  2-0  transition  an  additional 
contribution  to  arises.  The  latter  contribution,  however,  may  be 

neglected  in  virtually  all  cases  of  interest  because  it  originates  in  a 
polarization  with  spati  1  variation  other  than  that  of  the  inducing  field. 


Specifically,  P^  (z,  t),  the  velocity -averaged  polarization  at  frequency 

induced  by  a  pair  of  standing-wave  fields  E^(z,t)=Re  \  E°  cos  z  exp  i£2 
and  E9(z,  t)  =  Re  J  E?  cos  k„z  exp  ifi„t  t  is  of  the  form^, 


Pj  (?,t) 


Re  -i.  X  cos 

i  a 


k.  z  +  X,  cos  (2k0-  k. )  z  1 
lb  2  1 


ifi.  t. 
E°e  1 


where  X  and  X.  are  the  complex  susceptibilities  associated  with  the  response 
at  The  former  term  is  dae  to  the  interaction  of  pairs  of  travelling- 
wave  components  as  discussed  above,  and  may  be  obtained  in  our  treat¬ 
ment  from  the  sum  of  the  ensemble-averaged  polarizations  x  (v,  e  =  +  l) 
and  x  (v,  e  =  -  1)  lEq.  (C14)]  integrated  over  velocities.  The  latter  term, 
in  contrast,  is  due  to  the  coupling  of  several  travelling- wave  components 
and  only  becomes  large  when  and  are  tuned  close  to  their  respective 
atomic  center  frequencies.  We  are  thinking  of  a  polarization  induced  in  a 
sample  cell  placed  with  a  resonator.  The  net  emitted  power  may  be 
obtained  from  the  time  average  of  (z,t)E^(z,t)  integrated  over  the 
volume  of  the  sample  cell.  Consequently,  the  X  term  leads  to  a  contri- 

bution  to  the  emitted  power  identical  to  the  limiting  expression  of  Kq.  (58) 

2  2 

for  %  ^2  and  1^1  «  7  •  This  contribution  is  proportional  to  the  length 

of  the  sample  cell  and  is  independent  of  its  position  within  the,*  resonator. 


In  contrast,  the  contribution  due  to  the  term  vanishes  when  the  applied 
fields  are  detuned  several  natural  widths  from  their  atomic  center  fre¬ 
quencies.  Furthermore,  for  il  j  and  tuned  close  to  their  center  fre¬ 
quencies  the  latter  contribution  is  highly  sensitive  to  the  exact  location 
and  length  of  the  sample  within  the  resonator,  and  in  virtually  all  cases 
of  interest  is  either  negligible  for  all  locations  with  n  the  resonator,  or 
may  be  positioned  to  become  so.  Inspection  of  the  treatment  of  Ref.  2 
indicates  that  these  conclusions  are  also  applicable  for  cases  where 
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111.  APPLICATIONS  AND  DISCUSSION 


In  this  Pait  important  features  of  the  preceding  analysis  are 
illustrated  by  means  of  several  discussions  and  applications. 


_A  .  Frequency  Characteristics  of  the  Rest- Frame  Response 


It  will  be  recalled  that  the  Doppler-broadened  line  shape  was 
obtained  by  summing  the  ensemble -averaged  response  of  a  particu¬ 
lar  velocity  band  over  the  distribution  of  atomic  velocities.  This 
response  was  formulated  in  terms  of  and  [Eqs.  (33)  and  (34)) , 

I 

the  transition  rates  associated  with  the  emissior  of  photons  at  f2j. 

Let  us  now  inspect  these  quantities  in  the  atoms'  rest  frame,  prior 

to  Doppler  averaging.  To  avoid  unnecessary  detail,  consider  the 

2  2 

limit  of  these  expressions  for  weakly-saturating  E„  ( I  /3  I  «  y  ). 

(2  9) 

From  Eqs.  (35)  -  (38)  one  finds: 


J21(  n'r  fl!,)  =  21  a  I2!  3  !2Re-L  ! 


T0  i  i  i  * 


i  '  ■;  im  ? 

J01(  Slj,  «,)  =  2  |  or  r  - +  J,  ' 


3  2  2 

JQ1  =  2  1a!  &\ 3' 


Ll! 


01  , 


1  i  /  1 


(S  la) 


.  (61b) 


As  explained  earlier,  Jgj  is  the  2  -*  1  transition  rate,  due  to  double¬ 
quantum  exchanges  with  the  applied  fields;  the  first  term  of  JQ1  is 

the  0  -*  1  single-quantum  transition  rate  in  the  absence  of  E ^  and 

3 

the  second  term,  which  we  have  denoted  as  describes  the  nonlinear 
dependence  of  the  0  -♦  1  transition  rate  to  lowest  order  in  E^. 
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The  frequen  ',y  behavior  of  J9J  and  is  determined  by  the 

f  t 

quantities  L,  =  (fi  -  u, )  +  iy  ,  L„  =  (-ft.,  +■  w.,)  +  iy  and 

1  1  i  10  2  2  2  20 

t  i 

R  =  (ft.  -  w. )  -  (ft  -  w9)  -  iy  appearing  in  the  denominators  of 
1  1  ^  “  21 

the  various  terms  of  (60)  and  (61).  These  quantities  exhibit  distinct 
types  of  tunability:  *  and  *  may  be  associated  with  the  usual 

two  level  behavior ,  in  which  one  of  the  applied  fields  resonates  with 
a  pair  of  energy  levels;  R  \  on  the  other  hand,  involves  the  frequency 
of  both  fields  and  becomes  large  when  their  separation  approaches 
the  2-1  level  spacing.  Terms  involving  Lj  1  and  L9  *  together  are 
strongly  enhanced  when 

nrui  ■  ! 


simultaneously,  a  condition  equivalent  to  Eq .  (56),  as  discussed  in 
Section  IIB.  Terms  in  R  1  are  also  enhanced  under  those  conditions. 
Note,  however,  that  R  1  remains  resonant  for  the  less  stringent 
condition 


n2  ‘  “l 


u2  "  W1 


which  does  not  necessitate  resonant  coupling  of  the  two  level  type 

[Eq.  (62)].  Close  to  resonance  [Eq.  (62)]  J21  and  J^.  exhibit  different 

line-shape  characteristics.  Nevertheless,  in  this  region  they  are 

comparable  in  magnitude.  A  few  natural  line  widths  away  from 

resonance,  however,  J2j  begins  to  dominate,  becoming  larger  than 
8 

J01  by  a  factor  ~y/  !  Lj  ! .  In  fact,  in  this  limit  J<2  (  reduced  to  the 

(30) 

familiar  expression  for  the  Raman  transition  rate. 


■^21  ~*2  '  a 


>2,.o 


!  R !  2 1  Lj'2 


?,')  - 


and  becomes  negligible.  Thus,  double-quantum  transitions 
predominate  when  level  0  is  detuned  from  resonance.  These 
features  have  been  well  substantiated'  '  in  the  microwave  region, 
as  stated  above. 

In  considering  a  pair  of  applied  fields  resonating  with  a  broad 
atomic  velocity  distribution,  it  is  possible  for  atoms  within  a  par¬ 
ticular  velocity  band  to  satisfy  condition  (62);  consequently,  in  our 

3 

case  Jqj  and  contribute  comparably  to  the  velocity -averaged 
response. 

The  foregoing  discussion  is  directly  applicable  to  the  other 
energy  level  configurations  of  interest.  Note  that  in  extending 
these  considerations  to  the  cascade  case  [Fig.  2(c)] ,  the  "Raman" 

I  I 

condition  [Eq.  (63)]  should  be  replaced  by;  fij  +  n9  =  Uj  +  u>2. 

See  Appendix  C  for  further  details. 

So  much  for  the  frequency  response  of  a  particular  velocity 

band.  The  Doppler -broadened  response  is  obtained  by  summing 

this  quantity  over  the  entire  distribution  of  atomic  velocities.  As 

pointed  out  earlier,  in  the  present  case  these  averages  [Eq.  (46)] 

have  been  carried  out  in  the  fully  Doppler-broadened  limit  (y/ku«l), 

resulting  in  a  number  of  cancellations  which  greatly  simplify  the 

final  expressions.  Thus,  in  thir  limit  the  velocity  average  of  J2 

ii 

becomes  identical  with  the  negative  of  the  velocity  average  of  . 
This  may  be  easily  seen  ir.  the  special  case  of  weak  saturation  of 
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the  0-2  transition,  Eqs.  (60)  and  (61),  where  the  first,  term  of 

averaged  over  a  broad  distribution  of  velocities  vanishes;  the 

remainder  is  equal  to  -J  .  It  is  important  to  point  out  that  such 

2 1 

cancellations  do  not  occur  in  higher  orders  of  y/ku.  For  instance, 
the  complete  cancellation  of  y^  in  r^T  which  occurs  in  the  case  of 
(w./a^)  ~  1  (see  Section  IIIC,  below)  does  not  occur  in  the  next 
order  of  y/ku. 
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B.  Spontaneous  Emission  Line  Narrowing  Effect 

M  e  turn  next  to  the  manifestation  of  the  Laser-Induced  Doppler  I  ne 
Narrowing  Effect  appearing  in  the  profile  of  spontaneous  emission  from 

either  of  the  levels  of  a  Doppler-broadened  gas  laser  transition,  In  this 
case  the  interaction  of  a  standing- wave  laser  field  with  the  2-0  transition 
(Fig.  3  )  considerably  influences  the  0-1  profile  of  spontaneous  emission  as 
observed  along  the  laser  axis.  This  influence  arises  both  through  single¬ 
quantum  and  double- quantum  exchanges  with  the  radiation  fields  in  a 
manner  completely  analogous  to  the  induced  emission  case,  and  the  resulting 
line  shapes  are  identical.  In  fact,  the  0-1  spontaneous  emission  spectrum 
follows  the  spectrum  of  emission  induced  by  a  weak  monochromatic  probe 
field  tuned  through  the  resonance  when  the  population  of  level  1  is  ignored. 
Accordingly,  Eq.  (58)  gives  the  required  sponta- eous  emission  spectrum  when 
Nj  is  set  equal  to  0  andj^^i)  is  interpreted  as  the  usual  unperturbed  Doppler 
profile  (E.2=  0/  due  to  0  -  1  spontaneous  emission-  In  the  quantized  field 
calculation  of  Ref.  43  the  laser-induced  spontaneous  emission  profile  is  ob¬ 
tained  directly.  Those  results  are  in  exact  agreement  with  the  weak  satura¬ 
tion  limit  of  our  results  as  they  apply  to  spontaneous  emission,  given  below. 

To  amplify  the  foregoing  remarks  on  the  applicability  of  the  classical 
field  treatment  to  the  spontaneous  emission  case,  recall  the  discussions  of 
Section  IIA,  formulated  in  terms  of  the  response  (in  t he  atoms'  rest  farme)of  a 
particular  velocity  band  of  three  level  atoms  to  a  weak  monochromatic  probe 
field.  as  influenced  by  a  second  field,  strongly  interacting  with 

the  coupled  transition.  Suppose  that  the  monochromatic  probe  field  is  replaced 
by  a  wave  packet  distributed  in  frequency  about  over  a  narrow  interval 
dfi'j «  7  and  incident  in  a  given  direction  within  a  small  solid  angle 
dS:  Then  all  of  the  considerations  of  Section  IIA  follow  identically,  where  the 
field  intensity  c  E”  ^  8  sr  is  interpreted  as  an  average  quantity. 
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Since  df2^«  y,  the  response  at  is  independent  of  the  specific 
form  of  the  frequency  distribution  of  the  wave  packet.  Accordingly ,  let 
us  express  this  average  field  intensity  as 

—  E?2  =  I  d  £2'  (IS  ,  (65) 

O  7T  1  1 

where  I  represents  the  average  intensity  per  unit  frequency  interval  and 
solid  angle.  In  terms  of  F,  the  average  number  of  light  quanta  per  mode 
of  th(  (polarized)  radiation  field,  the  field  intensity  becomes 

ti«'  :3 

I  d  «  ^  cl S  =  P* — ~ ^dOjdS  .  (66) 

1  8  7T  '5  C  6 

The  transition  rates,  J.^,  representing  the  absorption  and  induced  emission 
processes  described  in  Section  IIA,  may  therefore  be  written  in  terms  ofT7  by  re¬ 
placing  E°  ^  in  Eqs.(33)  and  (34)  by  v  (ftoj  ^/ir^  c^)  dfl^  d  S  .  (Note  that 
enters  into  Eqs.(  33)  and  ( 34)  through  |a!^,  defined  as  |  q  j  q  E° /2li  |  )  To 

explicitly  emphasize  the  dependence  on  v,  let  us  write 

Jik=!r$ik"n,ids  <67> 

2  o  2  3  2 

with  $3ik=  (i/u  j  0 1  n^'/4h7r  c  )J  la!  .  The  effect  of  spontaneous 

emission  at  frequency  may  be  included  by  simply  replacing  v  by  j7  +  1 

(31) 

in  the  emission  rates,  JQ  ^  and  J^,  Consequently,  for  v  =  0,  the  l-*0  and 
1-*  2  absorption  rates  vanish,  and  the  0— ’  1  and  2-*  1  emission  rates  reduce 
to  dtt'^  d  S  and  $21^niclS’  respectively.  Using  Eqs.  (19a)  and  (19b) 
the  power  spontaneously  emitted  in  frequency  interval  du'^  and  solid  angle 


dS  is  seen  to  be 
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ti  «\  <n0?01+  n2&21)cU2l  tlS  •  (tH) 

Note  that  in  this  expression  the  n,  coefficient  has  vanished  because 
for  ~v  -  0  the  absorption  coefficients,  J?  and^.,,  also  vanish. 

vor  !  j'l  |  2  -*  0  Eq.  (68)  reduces  to  the  usual  expression  for  spontaneous  emission 

from  the  0  -  1  transition.  In  applying  this  formula  one  should  interpret  q 

(’  *  — ■*  a  3  — ’  . 

as  (R)  « j  •  ju  tp  ^  (R)  d‘  R  where  q  is  the  electric  dipole  operator,  t9  is  the* 

polari/.ation  of  the  laser  field,  and£j  is  the  polarization  of  a  spontaneously 
emitted  photon.  Thus,  the  angular  and  polarization  properties  of  the 
spontaneously  emitted  radiation  depend  on  the  characteristics  of  the  partic¬ 
ular  states  involved. 

Doppler  broadening  and  standing-wave  effects  may  be  included 
in  exactly  the  same  way  as  in  Sections  II 14  and  II  C,  respectively,  con¬ 
sidering  contributions  only  from  that  portion  of  the  spontaneous  radiation 
emitted  into  a  small  solid  angle  in  the  forward  (+7.)  direction.  Thus, 
as  stated  above,  the  spectrum  of  spontaneous  emission  from  the  0-  1 
transition  is  given  by  Eq.ta8)  with  N,  =  0  an interpreted  as  the 
usual  Doppler-broadened  spectrum  of  power  emitted  spontaneously  into 
dfl,  dS  with  given  polarization  in  the  absence  of  the  las  r  field.  Equation 
(58)  has  been  written  in  the  form  valid  for  both  folded  (a  +1)  and  cascade 
(or-  -  1)  cases.  Figs.  2(a)  and  2(c).  Suppose,  for  example,  the  laser  field  is 
tuned  to  a  frequency  below  its  center  frequency,  Wg :  Then  the  0-  1  spontane¬ 
ous  emission  change  signals  in  the  cascade  case,  which  result  primarily 
from  double-quantum  transitions  between  levels  2  and  1 ,  are  in  the  form  of 
resonant  increases,  with  the  1^.  resonance  above  Uj,  the  0-  1  center  f>-f  - 
quenev,  and  the  Tg  resonance  below  u>j  (Fig.  3).  In  the  folded  case, 
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on  the  other  hand,  the  change  signals,  which  primarily  result  from  single- 
quantum  transitions  between  levels  0  and  1,  are  in  the  form  of  resonant 
decreases,  mi  the  positions  of  F^  and  F^.  arc-  reversed  [Fig.  1(h)]. 

This  reversal  is  a  purely  geometrical  effect  arising  from  the  different 
way  in  wh  ch  (ho  fields  couple  to  the  respective  transitions  in  the  folded 
arid  cascade  cases.  (See  Appendices  G  and  D.  )  Similar  conclusions  were 
also  obtained  in  Ref.  4. 

Double -quantum  transitions  are  a  phenomenon  well  known  since 

/ on  \ 

the  early  days  of  the  quantum  theory.  Historically,  Eherenfest'1  'was 
ihe  first  to  suggest  a  model  capable  of  correctly  describing  sequential 
decays  among  levels.  This  model  was  motivated  by  the  correspondance 
principle  and  adiabatic  theorem  of  the  old  quantum  theory.  Results 
based  on  the  Eherenfest  model  were  later  found  to  be  in  agreement 

(33) 

with  the  quantum-electrodynamical  treatment  of  Weisskopf  andWigner. 

These  ideas  were  all  formulated  with  spontaneous  decay  in  mind.  How¬ 
ever,  in  the  presence  of  an  ntense  field  and  close  to  the  resonance 

>  i 

condition  idj  Uj  and  -  w.,,  an  additional  effect  becomes  important: 
he  intense  field  can  considerably  perturb  the  usual  single -quantum 
emission  rate  of  a  coupled  transition/* ***an  effect  comparable  in  mag¬ 
nitude  with  the  double -quantum  emission  rate.  Older  treatments  have 
not  been  concerned  with  the  former  effect;  line  shape  behavior  of  the  kind 

described  i.i  this  paper  involves  an  interplay  of  both  types  of  processes. 

(32  33) 

It  may  be  illuminating  to  outline  the  Ehercnh  st  picture.  ’  1  Consider 

the  cascade  system  of  l  ig.  2(c):  following  Eherenfest,  assume  that  when 
an  atom  >s  in  level  j  (energy  eigenvalue:  f.W.),  the  probability  that  its 
energy  lies  in  the  interval  between  f-x  and  Mx  -i  dx)  is  a  >rent/ian 
centered  aoout  ^W.  and  of  width 

i  1  t 


E^(x)dx  3 


(69) 


?r(Tj/2)dx 
(W.  -  x)“+(y./2 )" 

J  J 

(normalized  to  unity).  Consider  an  atom  initially  in  the  uoper  cascade 
state  (level  2)  which  underwent  decay  via  level  0  to  level  1,  accompanied 

»  I 

by  two  photons  successively  emitted  at  ft^  and  ft.,  (in  the  atom'  s  rest 
frame):  the  joint  probability  that  the  atom  had  energies  f  x,  f-.y ,  and  r  z  when 


it  was  in  levels  2,  0,  and  1,  respectively,  is  P  ,(x)P0(y)P  (z).  Evidently, 

»  I 

the  photons  emitted  in  that  event  were  of  frequencies  ft,-  x-y  and  ft^  y-z 

t  I 

Substituting  ft^  and  ft,^  into  this  probability  expression  and  integrating 
over  all  possible  values  of  y,  one  obtains  the  total  probability  of  an  atom 

1  I 

having  successively  errrtted  photons  at  i.j  and  ft 


fa 


(W2-ft2-y)2  -  (y2/2)2  (W0-y)2  +  <v0/2)2 


(Wj+ft'j  -y)2  v  (  y j /2)2 


(The  unimportant  numerical  factor  has  been  set  equal  to  unity.  )  Upon 
integration  (which  may  be  performed  as  in  Appendix  D)  Eq.  (70)  yields  a 
result  identical  in  frequency  dependence  to  the  lowest  order  expansion  of 
o  n  2  -*  i  emission  rate,  Eq.  (60)  or,  equivalently,  the  expression 
obtained  from  Eq.  (67).  (Equation  (60)  is  written  in  general  form.  See  Table  II) 

A  recent  discussion^ ,  formulated  on  the  basis  of  two -photon 
transitions  induced  by  a  weakly-saturating  laser  field,  has  analyzed 
the  frequency  profile  of  spontaneous  emission  arising  from  the  lower 
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laser  level  (Fig.  3).  This  discussion  is  presented  for  the  limiting  case 
•y9j-»  0.  It  is  easily  seer  that  in  this  limit  Eq.  (70)  involves  a  6-function, 

I  !  -2 

becoming;  6(^2  +  '  WjJILj1  ;  the  discussion  of  Ref.  5  is  equiva¬ 

lent  to  integrating  the  latter  expression  over  the  atomic  velocity  dis¬ 
tribution  for  t  -  +1  and  t  -  -1.  In  addition,  Ref.  5  states  the  line- 
shape  result  for  arbitrary  y^,  a  result  in  agreement  with  the  weak 
field  limit  (Q  «  i)  of  our  expression  (58)  as  it  applies  to  spontaneous 
emission  for  a  -  -1  and  Nq-0.  The  latter  result  utilizes  the  Eherenfest 

(34) 

model,  being  obtained  by  averaging  Kq.  (70)  over  velocities. 

S.»ch  an  analysis,  however,  ignores  the  important  role  played 
by  background  atoms,  Nq,  produced  in  level  0,  discussed  above.  This 
background  population  can  be  sizeable,  since  in  practice  the  populations 
of  upper  and  lower  laser  levels  often  differs  by  only  a  small  amount. 

Th?  extension  of  Eherenfest'  s  model  to  include  this  effect  does  not  follow 
in  an  obvious  way.  To  emphasize  the  significance  of  the  role  played  by 
background  atoms  in  level  0,  consider  a  cascade  system  in  which  only 
level  0  is  populated  (i.  e.  Nj  -  N9  -  0).  Then  in  the  rest  frame  of  an 

f 

atom,  an  applied  laser  field  at  H.  will  diminish  the  transition  rate  at 

I 

,  leading  to  two  holes  of  width  and  r^T  superimposed  upon  the 
velocity-averaged  emission  profile--an  effect  entirely  due  to  the  depen¬ 
dence  of  the  single-quantum  emission  rate  on  the  laser  field  intensity. 

On  the  other  hand,  the  analysis  of  Ref.  4  does  include  the  influence 
of  these  background  atoms,  although  that  treatment  does  not  make  explicit 
the  distinction  between  single -quantum  and  double- quantum  events. 


FIGURE  4:  Mode  Crossing  Effect,  Simplified  Experimemal  An 
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C.  Mode  Crossing  Experiments 

i  he  narrow  resonance,  ]  is  being  utilized  in  a  series  of  precision 

(12) 

measurements  of  hyperfine  structure  in  xi  non  and  G-fnetors  in  atomic 
(l  3'* 

oxygen'  ‘  *.  In  these  experiments  Wj  is  close  to  and  1^.  is  considerably 
less  than  rR-  Levels  1  and  2  are  a  pair  of  tunable  Zc*.  man  components 
optically  coupled  to  a  common  level,  level  0,  and  the  monochromatic  fields 
Lj  and  Lg  <*re  two  oscillating  laser  modes  of  fixed  frequency  separation  de¬ 
termined  by  the  cavity  length.  As  the  Zeeman  levels  are  magnetically  tuned, 
resonant  behavior  occurs  when  their  splitting  approaches  the  frequency  sep¬ 
aration  between  the  two  laser  modes.  We  shall  therefore  refer  to  this  effect*1  ^ 
as  "mode  crossing".  Important  aspects  of  mode  crossing  have  been  pre¬ 
sented  in  Ref.  2;  some  of  these  are  discussed  here,  for  completeness, 
together  with  additional  details. 

A  simple  version  of  the  mode  crossing  technique  is  shown  in  Fig.  4: 

The  output  of  a  Brewster-angle  gas  laser  oscillating  in  two  modes  is 
attenuated  (or  amplified)  as  it  passes  through  a  magnetically  tunable  external 
sample  cell.  The  respective  laser  modes  should  resonate  with  the  two 
coupled  Doppler-broadened  transitions  of  the  sample  gas.  For  example, 
in  cases  where  the  laser  levels  themselves  consist  of  closely  spaced  tunable 
components  the  sample  cell  may  be  a  discharge  tube  containing  the  same 
gas  as  the  laser  tube.  The  intensity  of  the  transmitted  laser  beam  is 
studied  as  a  function  of  the  level  separation.  At  low  field  intensities,  where 
saturation  eftects  are  negligible,  the  attenuation  changes  slowly  as  the  levels 
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arc  tuned,  appreciable  changes  occurring  only  when  the  frequency  sep¬ 
aration  of  the  Zeeman  components  varies  by  an  amount  comparable  to  the 
Doppler  width.  However,  when  one  of  the  applied  fields  r'.^i  saturates 
the  transition  with  which  it  resonates  (2  -  0),  the  attenuation  of  the  coupled 
transition  (1  -  0)  undergoes,  in  addition,  a  sharp  change  when  the  appropriate 
frequency  condition  is  met.  The  expression  for  this  behavior  is  readily 
obtained  from  the  «  =  + 1  case  of  Eq.(50)  Let  us  specialize  to  the  limit¬ 
ing  case,  relevant  to  mode-crossing  studies,  in  which 


i«2-  ui  *  “  <<:  "f 


where  u  and  y  characterize  the  magnitudes  of  the  most  probable  speeds 

(u-'s)  and  the  natural  widths  (y.'s),  respectively.  Assuming  the  velocity 
J  J 

distributions  are  all  Maxwellian  at  temperature  T,  the  power  emitted  or 


absorbed  at  in  a  small  length  of  the  sample  gas  is  given  by: 


-,w>f 

*  TV*  A’ 


i‘N0  7°C^  ^2l‘W21*"  o' ^71  +  72C^ 


in  whiciij^ftj )  was  introduced  in  h’q.  ( 57)  ^21  n2~ni’  W21  "  w?_  w  ,  end 
the  other  symbols  were  defined  in  Section  II B.  The  above  expression 
predicts  a  Lorcntzian  decrease  of  width  ®  as  *ne  Zeeman  :’Pl*u*nK< 

w 21 >  approaches  the  mode  separation,  of  weak  saturation 

(Q—  1).  (72)  reduces  to  an  expression  obtained  in  Ref.  2. 

There  are  two  noteworthy  features  of  this  behavior:  First,  that  the 
mode  crossing  frequency  condition  depends  only  upon  the  separation  be¬ 
tween  the  laser  modes,  being  insensitive  to  their  individual  frequencies; 


FIGURE  5:  Frequency  Stabilit; 
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and,  second,  that  the  width  of  the  mode  crossing  resonance  is  th<  sum  of 
the  widths  of  the  "crossing"  levels  (including  power  broadening  effects), 
and  is  essentially  independent  of  the  lower  level  width,  Yq-  These  two 
characteristics,  which  follow  from  the  approximation  of  closely- spaced 
crossing  levels,  (71),  have  far-reaching  implications  regarding  the  re¬ 
quired  frequency  stability  and  resolution  of  the  mode  crossing  technique 
as  utilized  in  precision  spectroscopic  measurements.  In  connection  with  fre¬ 
quency  stability,  note  that  although  the  frequencies  il  j  and  1 of  a  free- 
rurning  laser  may  wanner  over  a  sizeable  portion  of  the  broad  Doppler 
profile.  'fa  (flj ),  during  the  observation  time,  their  separation  will  remain 
virtually  fixed  (see  Fig.  5  ),  thereby  insuring  the  stability  of  the  observed  signal. 
In  regard  to  the  frequency  resolution,  note  that  in  many  important  gas  laser 
applications  the  widths  of  the  upper  levels  are  much  narrower  than  that  of 
the  lower  level;  the  complete  cancellation  of  the  latter  width  in  Eq.  (72) 
therefore  makes  for  a  relatively  narrow  signal.  This  prediction  is  strik¬ 
ingly  born  out  in  mode  crossing  studies  in  xenon^and  oxygen^  where 
the  observed  change  signals  are  not  only  much  narrower  than  the  Doppler 
width  but  are  also  more  than  an  order  of  magnitude  narrower  than  the 
known  radiative  widths  of  the  lower  level.  It  is  emphasized  that  these  two 
simplifying  features  are  characteristic  of  the  approximation  of  elosoly- 

/o  r.) 

spaced  crossing  levels,''  '  Eq.  (71),  and  do  not  occur  in  cases  where 
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It  should  he  noted  that  Eq.  (72)  holds  for  the  case  in  which 
E ^  is  weak  and  E9  may  fully  saturate  its  transition:  the  ease  where  both 
fields  are  intense  is  not  included.  N-  vc-rtheloss,  the  simple  behavior  of 
that  expression,  as  well  as  Eq.  (50),  its  more  general  counterpart,  is  re¬ 
vealing  of  the  effects  of  high  order  saturation.  It  should  also  be  noted 
that  when  both  fields  weakly  saturate  their  respective  transitions,  the 

saturated  response  up  to  "third-order"  may  be  recovered  from  (50). 

2 

This  may  be  achieved  by  expanding  (50)  to  lowest  order  in  1/31  ,  and 
noting  that  there  must  be  symmetn  between  a  and  ,3.  Thus,  aside  from 
the  usual  slowly  varying  third-order  background  term,  the  third-order 
response  at  £2^  is  just  given  by  (50)  with  Q=  1,  with  an  analogous  expres¬ 
sion  for  the  response  at  This  result  is  in  agreement  with  Ref.  2  . 

As  an  example,  Fig.  6  depicts  a  mode-crossing  signal  observed 

/ijp  v 

in  xenon  at  about  20  gauss.  The  experimental  arrangement  was 

similar  to  the  diagram  of  Kig.  4  ,  exc  ept  that  to  enhance  the  signal-to- 

noise  ratio  a  small  audio  frequency  component  was  superimposed  on  the 

slowly-varying  D.  C.  axial  magnetic  field.  The  detected  signal  was  fed 

into  a  phase-sensitive  amplifier  tuned  to  the  modulation  frequency,  and  the 

output  was  recorded  as  a  function  of  magnetic  field.  Accordingly,  the  curve 

obtained  in  Fig.  0  is  actually  the  derivative  of  the  mode  crossing  resonance 

(.37) 

of  Eq.  (72).  The  laser  source  was  a  3.37/u  xenon  laser  with  modes  separated 
by  about  50  MHz.  The  crossing  levels  were  pairs  oi'  Zeeman  components  of 
the  upper  laser  level  coupled  to  a  common  lower  level.  The  I  meter  mag¬ 
netically  tunable  sample  cell  contained  a  xenon  gas  discharge  at  a  pressure 
of  lip  Hg,  so  that  collision  broadening  was  negligible.  The  observed 
c -factor  is  0.929;  the  observed  width  of  1/2  MHz  is  consistent 
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with  estimates  of  the  upper  level  of  the  3.  37p  transition,  and 
about  30  times  narrower  than  that  of  the  lower  level. 

It  is  sometimes  convenient  to  place  the  sample  to  be  studied  within 
the  resonator  of  the  laser  source,  a  technique  useful,  for  example, 
when  the  laser  transition  itself  consists  oi'  closely-spaced  tunable  levels. 

As  disc  ’ssed  in  Section  IIC,  the  standing-wave  fields  within  the  laser  cavity 
produce  an  additional  change  signal  of  width  rR,  due  to  coupling  between 
oppositely  propagating  travelling-wave  components  of  E  and  E^,  symmet¬ 
rically  located  about  the  center  frequency  from  the  mode-crossing  reso- 
nance.(See  Fig.  7.  )  Applying  approximation  (7 1)  to  the  expression  for 
rR  ,  Eq.  (55),  one  finds  rR  =  Q  +  -7qQ.  which  can  be  considerably 

broader  than  the  mode -crossing  resonance  (of  width  r^-). 

The  condition  for  resonance  is  P. +  w.>-  In  contrast  to  the  fre¬ 
quency  condition  for  mode  crossing,  this  condition  is  highly  sensitive  to 
the  absolute  frequencies  of  the  individual  modes,  so  that  a  small  instability 
can  cause  the  modes  to  drift  away  from  resonance  (Fig.  7);  conse¬ 
quently,  observation  of  the  latter  signals  would  require  absolute  frequency 

(30) 

control  as  in  Lamb-dip  experiments. 

Finally,  we  note  that  in  mode-crossing  experiments  performed  within 
the  laser  cavity  effects  due  to  saturation  broadening  are  enhanced.  The 
overall  behavior  of  these  effects  is  expected  to  be  similar  to  Eq.  (72),  i.  e.  , 
to  broaden  the  resonance  by  the  saturated  widths  of  the  crossing  levels. 
Indications  of  such  broadening  has  been  observed  in  several  instances. 

For  example,  in  mode-crossing  experiments  performed  on  the  xenon  3.  37|i, 

(lO) 

transition  within  the  cavity,  the  widths  observed  are  generally  somewhat 
broader  than  that  of  Fig. 6. 
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I)  Spontaneous  Emission  Line  Narrowing  Experiments 

The  spontaneous  emission  line  narrowing  effect  has  been 

{ 1  o ) 

utilized  in  measurements  of  isotope  shifts  and  line  width  para¬ 
meters^  in  neon.  At  the  time  of  the  initial  observations,  a  de¬ 
tailed  picture  of  the  line  narrowing  effect  predicting  differing  widths 
lor  rx  and  i~B  was  not  at  hand.  The  width  difference  was  not  observed 
in  these  experiments  since  in  the  neon  transitions  studied,  and 
differ  by  only  a  small  amount:  its  observation  requires  good  laser 
stability  and  high  finesse  Fabry-Perot  analysis.  The  observation  of 
different  widths  in  spontaneous  emission  in  neon  has  recently  been 

(17) 

reported  by  Holt  . 

E.  Concluding  Remarks 

In  concluding,  we  would  like  to  make  additional  comments  con¬ 
cerning  several  aspects  of  the  Laser-Induced  Doppler  Line  Narrowing 
Effect. 

In  the  present  discussions  the  atoms  have  been  assumed  to  relax 
by  means  of  radiative  decay.  It  should  be  noted,  however,  that  Eq.  (41) 

and  its  subsequent  average  over  velocities,  Ec.  (50),  also  correctly 
deserve  relaxation  through  hard  collisions  when  the  are  interpreted 
as  phase-disrupting  hard  collision  rates.  In  fact,  in  the  treatment  of 
Kef.  10,  which  agrees  identically  with  the  limiting  case  of  (41)  of  the 
y  ■  all  equal,  relaxation  was  introduced  via  a  single  collision  rate.  The 
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lye  j 

detailed  features  of  the  line  shape  in  that  limit  have  been  fully  verified'- 
in  the  microwave  region,  where  Doppler  effect  is  negligible  and  the 
linewidths  are  entirely  due  to  collision  effects. 

On  the  other  hand,  in  situations  where  soft  collisions- -collisions 
in  which  phase  disruption  is  incomplete- -play  an  important  role,  additional 
details  not  considered  here  further  influence  the  line  shape.  Treatments 
such  as  Ref.  40  may  be  readily  extended  to  this  problem.  Inclusion  ol 
such  effects  is  highly  desirable  in  view  of  the  fact  that  the  pressure 
dependence  of  the  Laser-Induced  Doppler  Line  Narrowing  Effect  is 

observable  in  a  number  of  rather  different  experimental  contexts.  The 
classical  field  treatment  presented  here  seems  to  lend  itself  more  readily 

to  inclusion  of  collision  effects, as  opposed  to  quantized-field  treatments. 

The  line -narrowing  effect  presented  here  is  also  potentially 
applicable  to  high-resolution  microwave  and  optical  studies  of  inhomo¬ 
geneous  ly -broadened  solids.  In  gases  the  Doppler  effect  causes  the 
frequency  of  an  applied  travelling-wave*  field  10  appear  different,  for 
atoms  of  different  velocities.  In  solids,  on  the  other  hand,  the  fre¬ 
quency  of  an  incident  field  appears  the  same  to  a1!  of  the  atoms;  instead, 
inhomogeneity  of  the  crystaline  fields  may  result  in  a  broad  distribution 
of  atomic  center  frequencies.  The  procedure  for  including  inhomogen¬ 
eous  broadening  effects  is  similiar  io  the  Doppler -averaging  procedure 
employed  above.  However,  since  f/x^d  atoms  do  not  discriminate 
between  relative  propagation  directions  of  the  applied  fields,  line¬ 
narrowing  experiments  in  solids  will  produce  only  a  single  change 
signal,  the  analog  of  the  t  =  +1  change  signal  observed  in  a  Doppler- 
broadened  gas.  Another  important  distinction  is  that  the  simple  form 
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ot  the  relaxation  rates  introduced  into  the  Scht  bdinger  equation  is  no 
longer  valid,  since  non-radi&tive  relaxation  mechanisms  generally 
dominate  in  solids.  These  may  be  introduced  in  a  density  matrix 
formalism  as  in  Xnpendix  C  by  considering  appropriate  decay  constants 
for  the  diagonal  and  ofi-diagonal  density  matrix  elements.  The  inclusion 
of  h  effects  will  somewhat  modify  the  characteristics  of  the  observed 
change  signal. 
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APPKNDIX  A:  SVMMKTKV  UKTW  KKX  TRANSITION  PROP  ABILITIES 

Ft  is  state  I  in  Section  II  A  that  |c.  (t ;  t^(  k ) I  ,  the  transition  pro¬ 
bability  of  an  atom  produced  in  level  k  at  t0  to  level  j  at  a  later  time  t, 

2 

is  equal  to  I  c,.(t;  tQ;  j)!  ,  the  reverse  transition  probability. 

This  symmetry  may  be  demonstrated  by  direct  calculation.  Wo 
begin  by  restating  the  perturbation  technique  of  Section  II A  in  a  more 
formal  manner.  The  equatioi  s  of  motion  for  the  d.  (t,  tG,  k)  are  given  by 
(24).  To  obtain  a  solution  for  all  values  of  /land  jo-j  1,  the  d.  may  he 
expanded  in  powers  of  o\ 


dj  (1;  tt),  k)  =  dV  (t;  t0,  k)  +od,.(t;t0,k)  +  a  2d"  (t ;;  tQ,  k)  +  ...  (A  1 ) 


where,  as  before,  the  superscript  u  designates  the  parameters  of  the  un¬ 
coupled  system  (or  =  0).  We  require  f  solution  complete  to  Ofo),  con¬ 
sequent  only  the  two  leading  terms  of  (A  1 )  need  be  retained.  Inserting 
these  into  (24)  and  equating  coefficients  of  like  powers  of  a,  one  obtains 
the  following  sets  of  equations: 


"iU  l62{  u 

c'0  =  l[i(' 


dg  -  i/l^  e  “  d 


*  "  ^2  1  ,u 

0 


(!9 


dj  =  0 

V  i  l>  e  ‘  61  1  ,u  .  ,  1  62  1  ,1 
ci0  =l«  «<!  +  'm* 

d2  r  i3  c  “  dQ 


•i  .  "16 , t  ,u 

dj  5  *  1  do 


(A  2a) 
(A  2b) 
(A  2c) 
(A  3a) 
(A  3b) 


(A  3c) 
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Equations  (A  2a)  ami  (A 2b)  are  tin*  usual  equations  of  motion  a  damped, 
two- level  system.  (A  2c)  is  the  equation  of  motion  of  a  third,  decoupled  level. 
Equations  (A3)  describe  coupling  effects  to  lowest  order  in  a. 

From  (26)  and  (Ala),  the  boundary  conditions  are  found  to  be 


d^1  (t0:tc,  k>  =  d.  (t0;t„,  k) 
dj  (t0;  t(),  k)  =  0 

Thus,  for  an  atom  produced  in  either  level  0  or  level  2, 
(A  2)  and  (A  3)  reveals: 


(A  4) 


inspection  of 


dj  (t:  tn,  k)  =  0 


d.(t:t0,k)  =  0 


k  =  0,  2 


k  =  0,  2  ,  j  =  0,  2 


(A  5) 


Kq.  (28)  of  the  text  follows  immediately  fr<~>m  (A  3c). 

For  an  atom  produced  in  level  1,  inspection  of  (A  2)  reveals: 


dj  (t;t„,  1)  =  d.  (t0;t0,  1) 


(A  6 ) 


i6  jt 


Equations  (A  3a)  and  (A  3b)  may  be  eombim  d  to  yield: 

‘l0-i62‘,o+  ^!2(10  (62-  '  (A?) 

which  may  l)e  solved  subject  to  (A  4).  Then,  noting  (A  1 )  and  (A  6),  one 
obtains  <1q  (t :  *o>  1),  complete  to  0  (o): 

d„  {  sc  i6>'I\^.(ci'4T-,i<l-T) 

(A  8) 

*M+T  >q.T  y 

+  (q_e  -  q+e  )  V  . 


s[6\  V  61)  +  j 


The  symbols  employed  above  have  been  defined  in  Table  I.  Equation 


(A  8)  may  be  rearranged  to  yield: 

i n_T  i  n+T 

i6  t  q+  le  -  1]  q  [o  -  lh 

d0(t;to.  1)=|di(t0;t0,1)e  - - - ~  -  1 


-.(A  9) 


t  M.  9+  i 

Inserting  (A  9)  into  (A  3b)  and  integrating,  and  noting  (Al)  and  (A6),  one 


obtain? 


(&, -Mt  r 


in  T 


,(t;  t„,  1 )  =  dj  (tf):  t0,  1 )  e  “  |- 


-  1  e 


iq_T 


9+ 


—  l  .  (A  10) 


Finally,  utilizing  Eqs.  (15),  (16) ,  (26) ,  (39)  and(30),  the  symmetry  between  the 

probability  amplitudes  may  be  displayed; 

,  .  ,  .  -i  (Sl\  -  w.  )  t 

e ^  (tQ;  tf),  1 )  c  1  1 


cn(t;  t0.  1 )  =  Z-  - — -  — — , - —  c  (t:  t  C) 

0  o  CO^o'^o-0)  >(ul-wi)t0  1 


(Alla) 


e.,(t 


/4  1  1  \  ^  <)1  ”^1)1  )  1 

C1  to’ 1  e  -1  , 

;t0.l)=®-£  -2 -  - - T -  ci  (t:  t,,«  2) 


.,v< ,  <■(), 


ar|3  e9(tf  t0.  2)  (i-i(«2l-w2l)lo 


1 "  ’  lo’ 


(A  lib) 


Notice  that  from  these  probability  amplitudes  on  obtains 
(e0(t;to,  1 ) j2  -  |c1(t;l„,0)|2  , 

I  c20.:to,  l)j2  -  /Cj  (t :  t(),  2)|  2 


(A  12) 


as  anticipated  in  Eq.  (20). 


In  this  Appendix  the  relationship  ^ween  induced  polarisation 
an(1  transition  rates  is  examined.'"'  As  a  convenient  starting  point,  we 
reformulate  the  damped  equations  oi  motion.  Eq.  ( i  D.ln  •'  rms  ° 
density  matrix  with  elements 

,  ,  ,  .  (B 1 ) 

p..(t.to,k)e  <t  t  tf)>  k)  rj  (t:tn.k) 

,n  keeping  w  ith  previous  notation.  Pjj  ft;  tQ,  k)  denotes  the  value 
of  p.  a,  time  t  for  an  atom  produced  in  level  k  at  a  prior  time  t  and  the 
e s  are  the  probability  amplitudes  inroduced  in  Eq.  tub  The  familiar  equa 
lion  of  motion  for  p(t;t0,k)  follows  directly  from  (Bl).  (7)  and  (lit. 

,  i  (B  2) 

p  =  -i  <r.pl  -  j  '"'O' 

Here  [  1  are  commutator  brackets,  (  )  are  anti-eommutator  brackets, 

(B  3) 

V  V*i 

and  II  is  the  total  Hamiltonian  [Eq.  0)1 »  having  matrix  elements 


(B  2) 


(B  3) 


(B  4) 

15i.f*wiVvu  ’ 

with  V •  defined  by  (10).  The  expectation  value  of  an  operator,  such  as  the 

l.| 

induced  dipole  moment,  q  ,  is 


p =  tr  <p  p) 


<B  5) 


In  cases  such  as  the  present  one,  in  which  the  applied  fields  couple 
resonantly  lo  different  transitions,  the  equations  of  motion,  (B2I,  maybe 
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cast  in  another  form.  Inspection  of  the  equations  for  the  off-diagonal 
elements  p1Q  and  p2Q  reveals  that  in  the  present  ease  the  solutions  of 


interest  are  of  the  form 


^  ft  "  to<'<) 


in '  t 


m 


m  =  1 ,  2 


(116) 


wherf-  the  Am  are  slowly- varying  decaying  fund  ions  of  t  -  tQ.  Furth 


lermore. 


pjJ  arc 


inspection  of  the  equa.ions  for  the  diagonal  elements  shows  that  the 

(AO  \ 

likewise  slowly -varying  decaying  functions  of  t  -  t0.  iSee  the  discussion 
of  Appendix  C,  Section  2.  These  assertions  may  also  be  verified  directly 
by  forming  the  p..  (t;t0,k)  from  Eqs.  (27) .  ( 29).  ( 30)  and(All).  ]  Then,  from 
(B5),  the  dipole  moment  induced  at  frequen  y  S2  '  is 


m 


pm  *t:  to»  k)  "  p0mpm0  to*  ^  +  f'-  c- 


and 


ifi’  t 


Pm  ft;  tQ,  k)  =  i  «;T)  (wQm  pmf)-  c.  c.  )  +  (pnm  Am  e  m  c.  e.  )  . 


0  m  '  m 


(B  7) 


(B  8) 


Because  of  the  slow  variations  of  the  A^,  the  second  term  on  the  right  hand 
side  of  Eq.  (B8)  may  be  conveniently  neglected.  ^  Using  this  fact,  the 
diagonal  elements  of  (B2)  may  be  written  in  the  suggestive  form: 


Pll  (t  -  tf),  k)  +  Yj  pn  (t  -t0,  k)-  pj  (t;  t0,  k)E’  (/,  t), 


p99  (t-t0,  k)+ y9p  9  (t-tOJ  k)=  -_i_p  (t;t  k)E’(z,t)  , 
-  tifi;  c 


Pqo  ft  "  *o*  k)  +  Yq  Pqo  ft  “  lo»  - Pj  ft;  tn,  k)  E*  (7,  t) 


tin 


+  ^Li)2(t:t'i'k)E'(l'-t)  • 


(B  9a) 

(B  9b) 


(B  9c) 


II  is  convenient  at  this  point  to  introduce  sev.  ml  definitions  and 
identities.  I  he  initial  conditions  for  diagonal  and  off-diagonal  elements 
0l',)  lollmv  IVom  (lo)  and  (7)  and  may  be  written- 


PM  (t°:  {°'  k)  =  5ik6k, 


(B  10) 


Next,  note  that 


i-  r  Pjj(t'tO'k)Cjto  =  -  (  TT-P;;<t-t0,  k)dt 


J  at0  j.i  °’ 


Pi,  ^  -  to  =  o,  k)=-  6  , 

J  J  kJ 


(Bll) 


by  vn  tue  of  the  boundary  conditions.  Furthermore,  as  in  Eq.  (22)  we  in 

ti  oduce  (t),  the  ensemble-averaged  polarization  at  for  atoms  pro¬ 
duced  in  level  k: 


k  r1 

Pm(t)  =  nkTkj  pmft;t0,k)<lt0 


(B  12) 


Taking  the  time  derivative  of  (012)  and  noting  that  (mi)  implies  „  (t  ■ , 

'  m'  o- 

one  obtains: 


k)=  0, 


Pk  (t)  = 
m 


i>  (l:i„,k)dtf 


(B  13) 


I.et  us  inspect  Eq.  (Bite)  as  it  applies  to  atoms  produced  in  level 

21  k  =  2).  *  -uplying  both  sides  of  (B9a)  by  n^dt,,  and  integrating  from 
-  *»  to  t  using  (Bll)  and  (B13),  one  obtains 

n2  7 2  7  1  j  p  l  i  (t  -  lo»  2 )  d  tf)-  -  P“  (t. )  E'  (■/.,  t)/f»  H 
-  </"  * 

1  2  t  (B  1 4) 

"h  S2'  ^P1  *t)hl  (/,t)  \  +  (rapidly  varying  terms). 


'hlS  ut*ua,ion  tho  notation  t  ^  has  been  introduced  to  denote 


the  ti  me  average. 


-til  _ 


On  till'  right  hand  side-  of  (Bi  t)  \\  e  retain  only  the  i).  C.  term  --  tin-  remaining 
terms  art-  not  significant  and  wo  1  avc  systematically  ignored  contributions 
. (42)  !_>  ,,v  nrr.Cwi.vn  fra  r-  tltn  lr  — +  i  t  r-i*  Ci  i  t »  rtn  in- 


of  this  order. 


H  calling  ’lu-  expression  for  the  k  —  j  transition  i  ate  in¬ 


tv  jduced  in  Kq,  (18), 


Jkj^Yk7j)  k)dt° 


(B  15) 


(B14)  yields 


J  s.  J-.<pf(t)F.!(z,t)>)  . 
21  1  1  1 


(B  16) 


The  quantity  appearing  on  the  loft  hand  side  of  (B16)  was  formerly  obtained 
by  considering  the  probability  of  finding  an  atom  produced  in  level  2  at  tc  in 
level  1  at  a  later  time  t,  due  to  transitions  induced  by  the  applied  fields.  The 
right  hand  side  of  (B16)  is  the  rate  at  which  energy  is  emitted  at  frequency  in 
units  of  hflj  for  atoms  produced  in  level  2.  Accordingly,  -J91  represents 

the  rate  of  photon  emission  at  12 '  tier  atom  produced  in  level  2. 

& 

Next,  consider  Eq.  'B9c)  for  k  =  2.  A  procedure  similar  to 
the  above  results  in  the  expression 


V*o’-^r<*f'l>Ei<s’..tl-pTU  -P^OEju.t)  ;t  ( 

»  *2 

Combining  (B14)  and  (B17),  one  obtains; 
n9(.J9  +.J  )  =  — ! —  -  P9(t)E9(z,t)-  .  ( 

“  “  *  *  o'  -  -  t 

‘*2 

The  quantity  appearing  on  the  left  hand  s'de  of  (B18)  was  formerly- 
obtained  by  considering  the  probability  of  finding  an  atome  produced  in 
level  2  at  t  in  either  of  the  other  levels  at  a  subsequent  time,  due  to 
induced  transitions.  The  right  hand  side  of  (B13)  is  the  rate 


(B  17) 


(B  18) 


-  62  - 


at  which  energy  is  emitted  at  ft'  in  units  of  for  atoms  produced  in 
level  2.  Accordingly,  +  J2Q  represents  the  net  rate  of  phc.on  emis¬ 
sion  at  £*;  per  atom  produced  in  level  2.  Noting  that  J2J  is  associated 
with  emission  of  photons  at  V  [see  Eq.  (B16)j  and  also"  contributes  to 
photon  emission  at  ft,  ,  this  emission  rate  evidently  involves  exchanges 
ot  two  quanta  with  the  applied  fields.  This  identification  leaves  j  as 
tht  contribution  arising  from  single  quantum  exchanges  with  E^. 

Equations  (B16)  and  (B18)  were  derived  for  atoms  produced  in 
le%ei  2‘  Slmila^  manipulations  may  be  readily  extended  to  the  other 
possible  initial  conditions.  One  then  obtains  the  general  relations: 


mUk Jkm~  "  ^m(t)  Em(z*  l)>t  +  ^mV'k  6km*  m=l,2; 


V 


)  JL  i  "  % 

'  ki  -k 


(N'OIO  that  Eqs.  (B19)  have  beer,  derived  under  conditions  where  both 
fields  may  fuUy  saturate  their  respective  transitions.)  Combining 

these  equations  for  n.  I  and  using  («,.  one  finds  that  the  net  power 
emitted  at  ft.  is 


( B19a) 

( B19b) 


Iv(n|>=fia!  u,}  ,j.} 


1  "2  '’2r  no  'or  ni  (J 


1  WI2rJ10) 


(320) 


which  is  identical  to  Eqs.  (19)  of  the  text. 


APPENDIX  C: 
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CALCULATION  OK  POLARIZATION  USING  ENSKMD LE¬ 
AVER  AG  ED  DENSITY  MATRIX  FORMALISM 

This  Appendix  rmlorives  tin-  results  of  Section  IIA  of  the  text  by 
means  of  the  ensemble-avoi  aged  density  matrix  formalism  and  extends  those 
results  to  the  other  coupling  configurations.  Section  C  1  establishes  the 
connection  between  the  time-dependent  wave  functions  and  the  ensemble- 
averaged  density  matrix  equations  of  motion.  In  Section  C2  these  equations 
are  solved  for  the  "inverted-Y"  configuration.  Fig.2(a),and  the  ensemble- 
averaged  polarization  obtained;  the  same  level  system  was  treated  in  the 
text  by  means  of  the  transition  rate  approach  where  an  equivalent  result 
was,  of  course,  obtained.  In  Section  C3  the  method  of  Section  C2  is  ex¬ 
tended  to  the  other  level  configurations. 


1 .  Equation s  of  Motion 


The  elements  of  the  ensemble-averaged  density  matrix  corresponding 


to  Eq.(Bl)  of  Appendix  B  are  defined  by 

t 

Pij(r,  t)-  2nkYk  \  Pjj  <**>  V  k) 


(Cl) 


( Note  that  as  we  have  defined  them  here,  the  elements  of  0  {t;t  ,  k) 
are  dimensionless  while  the  elements  of  p  (T,t)  have  the  dimen¬ 
sion  of  number  of  atoms. )  Using  the  relationship 

dt  3  Pij  (t:to’  k),1,o"  Pij(T  ;t  *k>+  Jpp  (t;l0»  k)  dtQ  ,  (C2) 

-  jC 

and  the  fact  that  in  the  moving  frame  H,  the  total  Hamiltonian,  is  independent 
°t  tQ,  (B2)  leads  to  the  familiar  differential  equation:^*4* 

p  ( r ,  t)  =  -  {(p°-  p).  r  }  -  I  (  H,  p]  .  (C  A) 


.ill'll 
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Note  that  in  <C2)  the  appropriate  boundary  conditions,  from  (12),  are 

p. .  (t  •  t  ,  k)  =  6  ..  6 
Hij  o  o  ik  kj 

so  that  in  eq.  (C3), 


(C  4) 


p?.  (r,  t)  =  n.  5. . 
M  ij  »  6iJ 


tr  [v  p (r, t)| 


(C  5) 


with  n-  the  steady-state  background  population  of  level  j.  Using  ef  and 

(Cl),  the  ensemble-averaged  polarization  is 


(C  6 ) 


2.  Ensemble- Averaged  Polarization:  ''inverted-V"  Configuration 

We  procede  to  calculate  the  ensemble-averaged  response  of  a 
particular  velocity  band  of  atoms  whose  level  structure  is  in  the  form  of 
an  "inverted- V"  configuration,  as  shown  in  Fig. 2(a).  From  Eq.(C6),  the 
polarization  induced  by  F'^(z,  t)  in  the  moving  frame  is 

I  ifiit  , 

2Re(n01  p1Q)-  He  'x  (v,c)  2  Aj  (z)c  (C7) 

which  defines  x  (v.  c  ).  the  complex  susceptibility  for  the  atoms  moving 
with  axial  velocity  component  v. 

In  the  steady  state  the  p..-  0  except  for  population  fluctuations  which 
are  entirely  negligible  when  the  Under  these  conditions  the 

density  matrix  equations  may  be  compactly  rewritten  in  the  following  form: 


i\q=  r?0+  — -  (-!^*(i*i)VoI]-=.c. 

10  10  fi  i  yf,  yi  yn  1U  U1 


1 20  "  *  20 


70  7i  7o 

—  I  { . +  (j.  +  _L) 

7q  ^0 
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P?Q  “  c*  r* 


(C  8a) 
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.01  H01  f, 


■  T  (m01  rl0  +  p02  p21) 


^20P2C  =  +“  <p20r20  +  p10p21 


)  , 


(C  8b) 


X21  P21 


“  (m20p01  “  P01  P20J 


In  the  latter  equations  E*  =  E|  +  is  the  net  impressed  field  as  seen  in  the  atom’s  rest 
frame;  we  have  also  introduced  the  operator 

j i  —  -  (W-- W.)  +  i  Y- -  ,  (CO) 

dt  1  1  '1 

and  have  written  r.  =  p..  -  p  y  and  r?.  -  n.  -  m. 


This  set  of  equations  may  be  solved  by  inserting  e'  in  complex  form. 
The  various  Fourier  components  of  E'  drive  the  off-diagonal  elements  p.^. 

Of  special  importance  are  the  coefficients  of  the  p.  on  the  left-hand  side  of 
Eqs.  (C  8b).  These  arc  associated  with  the  resonant  behavior  of  the  induced 
polarization.  The  only  important  frequency  components  of  the  p.  are  those 
which  can  reduce  the  magnitudes  of  their  coefficients  to  the  y—  for  particular 
values  of  12  ^  and  12^.  Tn  other  words,  a  nearly  exact  solution  of  Eq.  'C8) 
can  be  obtained  by  assuming  the  following  form  for  the  off-diagonal  density 
matrix  components: 


o0i  = 


-iG  t 


n  t 

=  \e  2 


(CIO) 


^  i  («2  ~  “]>  1 
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in  which  X,  A.  and  D  are  constants.  Inserting  (CIO)  into  Eqs.  (C8).  one 
obtains  a  set  of  five  simultaneous  linear  equations: 


"  1  {r10"  r10*  =  !  "  X  +  (—  +  —  )oA 

'  yo  70  y\ 


+  1  ^r20~  r20*  “ 


-  ~  A  +  (-L  +  _l)/3\ 

y0  T0  72 


c.  c. 


-  c.  c. 


■L!  A=*r10  +  J3D  . 
L2X  =  2*  r20+Q,D 
R  D  =  o\  -  0*  A  . 


(Cl  la) 
(Cllb) 

(C  11c) 
(C  1  Id) 

(Clle) 


The  symbols  employed  here  have  all  been  defined  in  Section  II  A. 


A  solution  to  this  set  of  equations  for  all  values  of  E^  and  for  weak 

Kl  U-C-Iai<<y>  can  b''  Ohtamwl  by  means  of  a  simple  perturbation  technique. 
In  the  absence  of  coupling  through  u  the  system  reduces  to  a  simple  two-level 
system,  and  the  off-diagonal  matrix  elements  and  „2]  vanish.  Setting 
a  ~  0>  (C  1 1)  yields  the  unperturbed  solutions: 

O  q*  t  * 

Xu  -  r20>  ]"2 
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ru  =  ro  _  9  ‘20  'h  !  “  r9n 
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In  (C12)  the  superscript  u  has  been  introduced  to  designate  the  parameters 
of  the  uncoupled  system.  The  presence  of  coupling  through  Ej  does  not 
affect  the  unperturbed  parameters  X  ,  r^and  r.^Q  to  lowest  order  in  O'.  Thus 
the  first  order  coefficients  A*  and  L>'  are  determined  by  the  tirst  and 
third  off-  diagonal  equations,  (C  11c)  and  (C  lie): 


-L,  V  ■  A?0  +  3D’  . 


*  T-.'  *\U  I* 

R  D  2  a  \  -  p  A 


(Cl  3) 


Equations  (C  1  3),  together  with  (C  12),  (CIO)  and  (C  7)  yield  the  complex  en¬ 
semble-averaged  susceptibility  complete  to  first  order  in  u\ 

2  oL<>-2?2aR 


X (v, c  )  = 


"0,A-(„  -n  ('“lO1  R,,„  ...  x-io'  -r„ 


■■  tu.  —  1 1  r\  f  — — 

Aj(7.)  1  0  ft  B 

» 

The  time-averaged  power  emitted  atS2j  is 


+  (n  ^  -  nQ  1 


fs 


3  I 


AB 


1  v  ]  =  9  12  Ai  !  l,Tl  ^  <v*‘  > 


(CIS) 


the  present  result  is  thus  in  full  agreement  with  Kq.  (39)  of  Section  II  A. 


(C 14) 


3.  Extension  to  Other  Level  Configurations 


The  developments  of  the  previous  section  may  be  adapted  to  the  other  level 
configurations  in  a  simple  and  straightforward  manner.  The  configurations 
of  interest  are  illustrated  in  Fig.  2.  In  terms  of  the  relative  position  of  the  common 
level,  level  0,  they  are:the  "V"  configuration,  level  0  lowest  iFig.  2(b))  ;the  cascade 
configuration  ,  level  0  between  levels  land  2  [Fig  2(c)]  ;and  the  "inverted -V"condguration, 
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level 0  highest  lFig.2(n)].  The  equations  of  motion,  (CP,),  differ  in  the  various  cases 
only  in  the  relative  signs  of  \\\-  \Y.  which  enter  into  theoC_,  pjq,  (C  9). 

These  sign  changes  lead  to  different  choices  for  the  resonant  contributions 
to  the  off-diagonal  density  matrix  elements,  (CIO),  which,  in  turn,  neces¬ 
sitate  other  changes.  The  modifications  of  Section  C2  necessary  for  its 
extension  to  the  other  level  configurations  are  summarized  in  Table  II. 

Inspection  of  the  final  results  (noting  Ino  definitions  of  A*  in  Table  II)  re¬ 
veals  that  the  susceptibility  for  the  cascade  and  "inverted- V"  configurations 
art  the  same,  except  that  in  the  latter  case  A*  is  replaced  by  -A*  where- 
ever  it  appears.  Furthermore,  the  susceptibility  for  the  "V"  configuration 
is  the  negative  of  that  for  the  "inverted--),  "  configuration,  a  result  expected 
intuitively. 
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TAiil  K  II.  F.xtonsion  to  Other  Level  Configurations 
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(36a),  (36b) 
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"V"  [Fig.  2(b)]  ” Inverted- V"  l Fig.  2(a)j  cascade  [Fig.  2(c)] 
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APPENDIX  13:  VELOCITY  INTEGRATION  IN  THE  DOPPLER  LIMIT 

This  Appendix  presents  in  detail  the  velocity  integration  procedure 
outlined  in  Section  1IB.  As  seen  in  the  text,  Eq.(41>,  the  expression  for 
tne  Doppler -shifted  response  of  a  particular  velocity  band  of  atoms  viewed 
in  the  laboratory  frame,  is  a  rather  complicated  function  of  velocity.  More¬ 
over,  the  convolution  of  such  an  expression  with  an  atomic  velocity  distri¬ 
bution  G.{v),  Eq.  (47),  depends,  in  general,  on  the  specific  form  of  the 
latter  cistribution.  In  the  fully  Doppler-broadened  limit  y/ku«-  l,  however, 
it  is  possible  to  perform  the  averages  without  recourse  to  the  details  of 
the  G.'  s.  Furthermore,  a  number  of  cancellations  ocv.  ar  in  this  limit, 
resulting  in  a  rather  simple  final  expression  for  the  velocity -averaged 
response. 

1.  Velocity  Integration 

As  shown  in  Section  3  of  Appendix  C,  the  expressions  for  the  ensemble- 
averaged  response  for  the  various  level  configurations  of  interest  are 
closely  related  to  one  another.  It  is  worthwhile  to  perform  the  velocity 
averages  in  a  general  form,  applicable  to  all  of  these.  Referring  to 
Table  II,  the  susceptibilities  for  the  "V"  and  "inverted- V"  level  configu¬ 
rations,  Figs.  2(b)  and  2(a),  respectively,  differ  only  in  an  overall 
sign;  accordingly,  the  corresponding  velocity  averages  may  be  performed 
identically.  Furthermore,  the  latter  expression  may  be  readily  extended 
to  the  cascade  configuration,  Fig.  2(c):  following  Eq.  (46),  the  expression 
lor  the  Doppler-broadened  power  spectrum  may  be  written: 
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2l.Hjl.rrim 


c )  dv: 


(1M) 


the  subscript  a  has  been  added  to  denote  the  "invertcd-V"  level  configu¬ 
ration  (a  =  +1),  treated  in  the  text,  and  the  cascade  configuration (cr  -  -1). 
As  in  Eq.  (40), 
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(D3) 

(D4) 


From  Table  II,  the  appropriate  generalization  of  Eqs.  (36)  is 

Lj  =  a',  +  iY10  ,  (D5a) 

L2(^)' _cr  ^2  +  i "^20  ’  !D5b) 

R(a)  =  (A*  -ctA2)  -  i721.  (D6c) 

To  perform  the  averages,^ (v,  *)  is  decomposed  into  its  partial 
fractions.  In  addition  to  n^N.Gj(v)lsee  Eq.  (47) ) ,  J^is  velocity  depen¬ 
dent  through  Lj ,  L2,  and  R,  Eqs.  (D5).  Explicitly, 

ct=+1:  "inverted-V"  (see  text) 
a-  -1:  cascade 
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with  A  .  = 
J 


in  which 


and 


Also, 
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(D6a) 
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S2 .  -  cj..  To  expedite  the  decomposition,  note  that 
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and  Q  is  the  factor  by  which  the  saturated  levels  of  the  0-2  transition  are  broadened: 
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(OH) 


Inserting  Eqs.  (139)  and  (D12)  into  (i32),  the  partial  fraction  decomposition 
leads  to  expressions  of  the  form  of  Eq.  (48)  of  the  text: 

H 

•  I  I 

in  which u  ,  o  ,  y  and  y  are  real  and  y  and  y  >  Q.  As  stated  in  Eqs,  (49) 
of  the  text,  in  the  limit  y/ku  1, 
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A  proof  of  Eqs.  (D16)  will  be  found  at  the  close  of  this  Appendix.  Thus, 
in  the  velocity -broadened  limit,  the  question  of  whether  or  not  a  particular 
partial-fraction  term  of  the  typo  (III 5)  contributes  to  (Dl)  is  entirely 
determined  by  the  signs  of  the  imaginary  portions  oT  its  corresponding 
factors.  Many  of  these  are,  in  turn,  controlled  by  the  sign  of  K  (  Eq.  (D7)|. 
We  shall  confine  our  attention  to  the  important  case  in  which  K  2:0,  always 
valid  if  either  k  ^  >kv  or  at  =  -  1 .  The  complementary  case  in  which  K  '  0 
can  be  obtained  similarly. 

It  follows  from  their  definitions  that  when  K  0,  0. 

Using  these  facts,  the  velocity -averaged  partial  fraction  components  of  I 

(7 

may  be  combined  tc  yield: 
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n 

Also,  in  evaluating  the  G-'s  in(DlT)the  approximation  that  j0i  /y9  «  ku 
has  been  made. 

The  first  term  of(Dl?)will  be  recognised  as  the  power  emitted  by 

an  unsaturated  Doppler-broadened  two  'eve!  system.  The  second  term  may 

be  simplified  bv  utilizing  the  identitv(D9)with-<x  (k.  k9)m  substituted  for 

t  -  <T€ 

kj\  .  Then  usirsg(  D4)  and  (D6),  (D18)reduees  to: 
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which  may  be  further  simplified  by  noting  that,  inserting  (DM)  the  numerator 
is  a  factor  of  the  denomenator.  Replacing  k,  / k,  by  w,/m9,  slight  re- 
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a  remarkablv  simple  result! 


For  cr  -  +1  this  expression,  of  course. 


t  .* 


reduces  to  Eq.  (50)  of  Section  IIB. 


D2.  Proof  of  Equations  (DI6) 

Consider  the  integral 

Z  k  G .  (v)dv 

Zdw,  ku,  *,)  r  \ - 1 


(u  +  iy)  +  kv 


(D23) 


in  whichw  and  y  are  real,  y>  C,  and  G^(v)  is  a  slowly -varying,  even  function 
of  v,  of  width  characterized  by  u.  In  the  limit  7 /ku  <  '  1, 

(u  +  iy  +  kv)  1  ~»P(~  -t-  kv)  1  -  i  6(-^-  -!■  v), 

where  P  denotes  principle  value,  and  (D23)  becomes 
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For  compactness,  we  shall  write 

Zj(w,  ku)  =  Zj(w,ku)  -  iZj  (u,  ku) 

,  •» 

with  Z-  and  Z.  the  real  and  imaginary  parts  of  (D25),  respectively. 

i-  -f- 

Consider  now  the  function H^”(v),  Eq.  (D16):  Decomposing  IP- (v)  into 
partial  fractions  and  then  using  (D23),  (D25),  and  (D26),  one  finds  that  in  the 
Doppler  limit. 
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The  only  significant  contribution  to  (D27)  arises  from  the  immediate  vicinity  of 
!  1 

o  -  (k/k  )w  =  0.  Under  the  stated  assumptions  the  Z.  vary  slowly  over  this 

^  (45) 

small  region:  consequently,  it  is  permissible  to  evaluate  the  numerator  ’  at 

,  •! 

u  =  -  (k/k*  )u'  .  Making  use  of  the  fact  that  Z^  and  Z^  are  odd  and  even 
functions  ofw,  respectively,  this  substitution  leads  directly  to  Eqs.  (D16). 
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(2  ku) 


X  ((!h12!',/Y1Y2><VNl,El3ll  +  1'>2lT12‘i(“2!‘''l)l  ' 

*  (l»12i4/Y1)F,E22l(Y12/Y2)W2-N1)(Y12-i(“2r,'B)l  ‘ 
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4(l1i12P,1l2/2Y,)<N34N2-2Nl>E,3|Y-l<UB-h1))-1 

+  (|p12M31i2,2Y,)(S3-Nj)E,E22lY-i<v>B-fB)l  ‘ 


(33a) 

(33b) 

(33c) 

(33d) 
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-1 


P;=-[i’rl/2/(4ku)llMl2M;i]|2(N2-Ni)E1E 22lY,2- ifr>21- 


1y32-  i(w32-A)1 


-1 


(33e) 


(33f) 


It  should  be  noted  that  this  expression  only  holds  when  the  two  Doppler 
gain  profiles  closely  overlap  and  the  laser  frequencies  are  fairly  close 
to  their  respective  center  frequencies.  More  generally,  each  Lorentzian 
term  is  multiplied  by  a  slowly -varying  Gaussian.  For  example,  term 
(33c)  vanishes  when  there  is  no  overlap.  A  complete  Erratum  for  Ref.  2 
is  being  published  concurrently  in  Phys  Rev. 

2f)-  In  obtaining  these  results,  the  first  term  of  Eq.  (34)  must  be  expanded 

to  first  order  in  Ij3|  .  In  the  remainder  of  (34).  which  is  identical  to 

(33),  B~  may  be  replaced  bv  (RL*)'1,  and  A'1  bv  I  L„  !  "2  =  (i/9  -y  ) 

I  2  ' '  r20 

d1  l2)  +  ^  L2  ^  ’  aiul  sPeciftc  use  is  made  of  the  fact  that  L *-r  L*  = 

1  2 

R  “i Yq* 

30.  Note  that  when  Ej  and  Eg  are  detuned  from  their  respective  center 

frequencies  by  sizable  fractions --a  situation  which  does  not  concern 

u»  here -- additional  nonresonant  terms  which  were  dropped  from  the 

Schrodinger  equation  (24)  may  no  longer  by  ignorcfi.  For  example, 

we  have  ignored  terms  with  frequency  dependence-  of  the  type  (ft'  -  u,  )* 1 . 

1  2 
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44.  W.  E.  Lamb,  Jr.,  andT.M.  Sanders,  Jr.,  Phys.  Rev.  119,  1901  (1960). 
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FIGURE  CAPTIONS 

Fig.  1.  Spectrum  of  lino  narrowing  induced  in  the  Doppler  gain  (or  atten¬ 
uation)  profiles  of  (a)  laser  and  (b)  coupled  transitions.  In  the  energy 
level  configuration  depicted  the  effect  takes  the  form  of  depletions  in  gain 
of  the  coupled  transition  over  two  narrow  frequency  intervals.  As  an 
example,  for  a  weakly  saturating  laser  field  and  for  closely  spaced  lower 
levels  of  equal  lifetimes,  the  broad  change  signal  (width:  r^)  of  the  coupled 
transition  is  twice  as  wide  as  the  laser  change  signals  (width; r^),  while 
the  narrow  change  signal  (width:  r^)  is  reduced  by  twice  the  width  of  the 
upper  level  (level  0)  and  is  independent  of  it.  In  the  cascade  configuration 
(not  shown),  level  2  above  level  0,  change  signals  of  similar  widths 
selectively  enhance  the  background  profile  of  the  coupled  transition.  Also, 
in  the  latter  case  the  positions  of  the  two  change  signals  are  interchanged. 

Fig.  2.  Energy  level  configurations.  Frequencies  Uj  and  fig  refer  to  the 
applied  fields  and  and  Wg  to  the  level  separations.  In  each  case  level 
0  is  common  to  both  1-0  and  2-0  transitions. 

Fig.  3 .  Spontaneous  line  narrowing  effect:  simplified  experimental  arrange¬ 
ment,  energy  level  diagram,  and  spontaneous  emission  line1  shape  as  in¬ 
fluenced  bv  laser  field.  In  the  latter  diagram  the  laser  is  assumed  to  be 
oscillating  somewhat  below  its  atomic  center  frequency  (u-  n_>  y  +i.  y  ). 

fa  faU  2  ^ 


See  text  for  details. 


~  *’  M°d0  ”rosslnS  effect,  simplified  experimcn  al  arrangement  and 
energy  level  diagram. 


Frequency  stability  of  mode  crossing  resonance.  The  frequency 
separation  of  the  applied  fields.  !i21,  remains  stable  against  frequency 
drift  caused  by  length  changes  of  the  laser  resonator.  Thus  the  resonant 
condition  n21  =  ui21  is  maintained,  regardless  of  drift  of  the  laser  fre¬ 
quencies.  Note  that  in  an  actual  experiment  the  Doppler  primes  may  overlap, 

^°de  crossing  change  signal  observed  in  Xc  at  3.  37p.  The  sample 
cell  is  external  to  the  laser  resonator  (Fig.  4). 

Fig^.  Frequency  instability  of  the  additional  standing-wave  resonance  which 
occurs  when  the  sample  cell  is  placed  within  the  laser  resonator.  The  condition 

for  the  latter  resonance,  0,^.  is  destroyed  by  slight  iluctuations 

in  resonator  length.  (Compare  Fig.  5. ) 
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Photoeloctnc  Potential  from  Photoreceptor 
Cells  in  Ventral  Eye  of  Limulus 

J.  E.  Brown,  J.  R.  Murray,  and  T.  G.  Smith 
ABSTRACT 

Intense  colored  light  from  a  CW  gas  laser  evokes  a  photoelectric 
potential  (PEP)  in  the  photoreceptor  cells  of  the  ventral  -ye  of  Limulus, 
This  PEP  has  two  components,  both  of  which  have  the  action  spectrum  of 
a  530  nm  rhodopsin.  The  evidence  is  consistent  with  the  hypothesis  that 
the  PEP  arises  directly  from  the  orderly  array  of  rhodopsin  molecules 
which  are  an  integral  part  of  the  photoreceptor  cell  membrane. 
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In  n  previous  paper  v.o  reported  In:*. I  when  the  lateral  cyt 
I.imulus  polypheuu  .  is  si  '.eu  1  ;\tui  with  very  uiUt-c  lights, 
a  short-latency  response  can  be  record  d  with  intracellular 
electrodes  X rum  the  eye’s  photoreceptor  cel  If.  (the  retinular 
cells)  (1).  This  electrical  response,  which  we  called  a  photo¬ 
electric  potent' al  (or  PEP)  .  was  found  to  have  a  number  oi  chnrnc- 
teristis  similar  to  the-  so-called  early  receptor  potential  (EKP) 
studied  in  vertebrate  eye  with  extracellular  electrode  .  (11-15)  . 

In  addition,  the  observations  that  the  polarity  of  the  pup  re¬ 
versed  across  the  phot ore \  lor  cell  i  »n-  and  th  amp]  i  tilde 

and  the  polarity  of  the  PEP  could  be  altered  by  chr.ncin"  membrane 
potential  indicated  that  the  {jonewator  of  the  FEP  was  intimately 
associated  with  the  photoreceptor  cell  membrane  (1).  However , 


there  remained  the  questions  of  whether  the  1, inulus  PEP,  like 
the  T  IP,  had  two  components,  and  ii  the  .spectral  sensitivity  of 
these  components  resembled  the  absorption  spectrum  of  Limn  1  us 
rhodopsir.  In  this  paper  we  present  evid  nco  which  answers  these 
questions  affirmatively. 

The  present  experiments  were  done  on  the  ventral  eye  of 
Limulus.  The  photoreceptor  cells  of  the  ventral  eye  have  many 
anatomical  properties,  membrane  characteristics  and  responses  to 
li^ht  which  arc  essential? v  the  same  as  those,  found  in  the  retic¬ 
ular  cell  of  the  lateral  eye  (17,19).  The  methods  used  in  thes> 
experiments  were  the  sane  as  employed  previously  (1),  with  two 
additions.  Here,  the  action  spectra  of  the  PEP  were  determined 


by  usin';  a  CV  pas  laser  (arp;on  or  kry contracted  and 


developed  by  one  of  us  (J.R-M-,  16).  With  this  laser ,  five  line 
of  sufficiently  ini*nse  monochromatic  light  within  the  spectral 
region  of  interest,  vt -re  available.  The  action  spectrum  of  the 


generator  potent  i  a  1  was  dote  mined  bv  uJ'.ing  C 

A- 

narrow-band  inter Lcr  nee  filters  and  neutral 
When  a  pulse  of  white  light  of  moderate 


Xenon  a>‘e  source 
density  filters, 
intensity  (ir.  the. 


range  of  mi  ero  .=; 


.)  is  focused  onto  a  ventral  eye  photoreceptor 


cell, 


the  membrane  potential 


undergoes  a  transient  depolariza¬ 


tion  vill*  the  c’>.  motorist ic  waveform  of  the  generator  potential 
(rig-1 ,A,  Vj)  and  with  a  latency  of  several  milliseconds  (Fig . 1A , 
cf.\'2  and  I.).  When,  however,  the  light  intensity  is  increased 
several  hundred-fold,  another  response  appears  whose?  onset  is 
coincident  with  the  light  pulse  (Fig. ID,  cf.V2  and  L) .  This 
response  is  the  PEP. 

From  Fig.  1,B,V2>  it  can  be  seen  that  the  PEP  has  two 
components.  The  first  is  a  depolarizing  potential  change  which 
is  coincident  with  the  lignt  pulse.  The  second  component  is  a 
hyporpolarizing  potential  change  which  outlasts  the  light  pulse. 


Monochromatic  light  from  the  laser  produced  the  same  PEP’s 


as  those  evoked  by  white  light.  The  action  spectra  of  both 
co  -portents  of  the  PEP  fit  reasonably  veil  with  the  absorption 


spectre.?  of  a  rbodopsin  (predicted  by  Dart  nail's  nomogram)  with 
the  absorption  max  rum  taken  at  530  nu  (20;  Fig.  2).  Previous; 
studies  have  shown  by  microrpect rophoto: entry  (18)  and  ifcith 
clectrophysiologicrl  techniques  (IS)  that  Jic  absorption  spectrum 
and  generator  potential  action  spectre  :,  respectively,  peak  near 
53C  n  We  cent  ir:;  that  the  ventral  eye  generator  potential  action 


spectrum  peaks  r  ar  530  rn  (Fig-  2). 


In  conclusion,  wo  have  found  that  a  PEP  can  bo  recorded 


with  intracellular  electrodes  in  the  photoreceptor  tills  of  the 
Li  mi  l  u.  ventral  eye.  This  PEP  is  similar,  in  t..o  more  respect. ••  , 
to  the  LHP  recorded  ex trace 1  1  ulerly  l ro; i  vertebrate  eyes,  viz.  , 
there  are  two  to;  poneats  to  t i>>  response-  and  both  components 
have  action  spectra  vhich  are  the  sane  as  both  the  absorption  spot 
truia  of  rhodopsin  and  the  action  spectrum  of  the  generator 
potential. 

We  and  other::  (2-15)  interpret  the  presently  ava.il able  evi¬ 
dence  to  indicate  that  the  PEP  and  VKV  ropre:  o  t  an  early  and 
perhaps  direct  eleetr  :.ic  change  in  rhodopsin  molecules.  Recent 
evidence  suggests  that  the  generation  o>  the  E!*P  is  dependent  on 


preservation  of  an  ord  red  array  of  rhoclopsin  molecules  (15). 
Moreover,  as  will  be  reported  1:  tor  (17),  the  properties  of  the 


ventral  eye  PEP,  like  those  of  the  lateral  eye  PEP  (1),  indicate 
that  the  PEP  generator  is  an  integral  constituent  of  the  photo¬ 
receptor  membrane.  Hence,  the  generators  of  the  LHP  and  PEP  may 
well  be  an  ordered  array  of  rhodopsin  molecules  lying  in  the  cell 


membranes  of  photoreceptors. 

Whether  the  PEP  and  EP.P  are  in  the  direct  chain  of  events 
leading  to  the  permeability  changes  associated  with  the  generator 
potential  and  ELI,  respectively ,  and  if  so,  how  they  arc  ccupltd 
to  these  permeability  changes,  remain  open  questions. 
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FIGURE  LEGENDS 


Fig-  1. 

Intracellular  responses  of  Limn  lus  ventral  eye  to  light  . 

In  both  A  and  B ,  Vj  is  low  gain,  D.  C.  -coupled,  slow  sweep- 
speed  recording;  V2  is  high  gain,  A.  C.  -coupled,  last  sweep- 
speed  recording;  and  L  is  light  monitor  synchronised  with  V, 

4U 

at  the  same  spcect*  In  i  t  .io  inrtral  pulse  is  a  voltage.-* 

time  calibration  i  rl.e.-  o>  1  mV  and  N  msec.  The  gain  of  the 
Vj  trace  is  20  tim's  less. ;  the  sweep  speed  5  tim.es  slower.  In 
A,  the  light  intensity  was  in  the  range  of  ±  watts.  In  B, 
the  light  intensity  was  in  the  range  of  m  watts. 


Action  spectra  of  toe  PEP  and  generator  potential.  The  relative 
sensitivity  (ordinate)  against  wavelength  in  nanometers  (na; 
abscissa).  Solid  line  taken  from  Dartnail  nomogram  for  530 
nra-pignent .  Open  circles  show  the  action  spectrum  of 
generator  potential.  Squares  and  triangles  indicate  the 
action  spectra  of  the  depolarizing  and  hyperpolarizing  com¬ 
ponents,  respectively,  of  the  PEP. 
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Motional  Narrowing  in  Hydrogen  Haman  Scattering 
J.  R.  Murray  and  A.  Javan 

This  report  describes  the  results  of  a  high  resolution  study  of  the 
motional  narrowing  effect  ^the  decrease  in  the  normal  Doppler  linewidth 
when  the  mean  free  path  becomes  comparable  to  the  wavelength  of  the 
emitted  radiation  (J_  -  5)^  .a  spontaneous  Raman  scattering  in  the  (^(1)  line 
of  the  1  -0  vibrational  band  of  the  hydrogen  molecule.  An  argon  ion  laser 
is  used  as  the  exciting  source,  and  the  linewidth  analysis  is  done  with  a 
Fabry-Perot  interferometer. 

We  have  reported  qualitative  preliminary  observations  of  this  effect 
previously  (6).  Motional  narrowing  was  first  seen  in  hydrogen  by  Rank  and 
Wiggins  in  infrared  absorption  spectra  (7).  Narrowing  has  also  been  seen 
in  stimulated  Raman  scattering  in  hydrogen  (8).  Cooper,  May,  et  al.  have 
obtained  similar  high  resolution  results  for  the  rotational  Raman  lines  of 
hydrogen.  (9) 

The  argon  ion  laser  operates  at  4880A  and  is  a  hollow  cathode 
design  with  segmented  graphite  discharge  tube  (10).  The  laser  is  operated 
in  a  four  mirror  Michelson  type  cavity  (II)  to  reduce  the  number  of  simul¬ 
taneously  oscillating  modes,  and  has  a  linewidth  of  about  0.  0085  cm  . 

(The  linewidth  under  similar  conditions  in  a  two  mirror  cavity  is  about 
0.  1  cm  * . )  The  center  of  this  line  is  locked  to  an  external  Fabry-Perot 
reference  to  reduce  frequency  drifts,  which  otherwise  are  a  severe  limit¬ 
ation.  The  output  power  is  about  50  milliwatts. 

The  laser  beam  enters  a  scattering  cell  which  is  a  quartz  capillary 
of  0.  5  mm.  inner  diameter  and  30  cm  length.  Reflection  at  grazing 
incidence  from  the  walls  of  this  capillary  propagates  the  light  scattered  in 
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a  small  cone  in  the  forward  (parallel  to  the  laser  beam)  and  backward  (anti¬ 
parallel)  directions  to  the  ends  of  the  cell,  where  it  is  collimated  by  a  lens 

4> 

_  1 

and  enters  a  pressure  scarined  Fabry -Perot  interferometer  of  0.  79  cm 

interorder  spacing.  The  output  of  the  Fabry-Perot  is  detected  by  an  ITT 

FW-130  photomultiplier  with  pulse  counting  electronics  and  displayed  on 

a  chart  recorder.  Recordings  are  made  at  a  temperature  of  25°C.  over  a 

density  range  of  2  to  100  Amagat.  Sample  recordings  exhibiting  the  line 

narrowing  effect  are  shown  in  Figure  1.  The  dependence  of  the  Q(l)  line- 

width  on  density  is  shown  in  Figure  2.  The  ristrumental  linewidi.  0.031  cm  *) 

has  been  subtracted  from  each  point. 1 

At  high  densities,  ignoring  pressure  broadening  effects,  all  models 

of  the  motional  narrowing  predict  that  the  normal  Gaussian  Doppler  profile 

D 

will  become  a  Lorentzian  lineshape  with  full  width  at  half  maximum  _ _ 

tr  d 

where  DQ  is  the  self-diffusion  constant,  k  is  the  wave  vector  of  the 

incident  photon  minus  the  wave  vector  of  the  emitted  photon,  and  d  is  the  density 

in  Amagat  units  (density/density  at  1  atm.  pressure  and  0°C).  If  pressure 

broadening  is  present  and  is  not  correlated  with  the  motional  narrowing 

D  k^ 

{£-5),  the  resulting  lineshape  will  be  a  Lorentzian  of.  width  +  ad, 

where  a  is  the  pressure  broadening  coefficient,  and  is  determined  from 
the  forward  scattering  data  to  be  1.  68  ±  .  07  x  10  3  cm  V^magat,  which  is 
in  fair  agreement  with  the  calculations  of  Van  Kranedonk  (12)  and  the 
measurement  at  higher  densities  by  May,  et  al.  (13). 

The  solid  line  for  the  backward  scattering  in  Figure  2  gives  a  fit  to 
this  simple  diffusion  model.  The  dotted  iine  gives  the  behavior  of  a  "hard 
collision"  model  0,4)  in  the  region  where  it  diverges  from  the  simple  diffusion 
model.  This  model  assumed  that  the  velocity  of  a  molecule  after  a  collision 
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has  a  Maxwell -Boltzmann  distribution  and  is  unrelated  to  its  velocity  before 

the  collision.  The  diffusion  constant  D0  has  been  chosen  for  a  best  fit  to  this 

model,  and  is  1.  35  ±  .  10  cm  Amagat/sec.  A  gas  kinetic  measurement  of 

Dq  (1_4),  expressed  in  our  units,  and  corrected  for  temperature,  gives 

D0  =  1.  36 1  ±  .  004  cm  Amagat/sec.  Asa  measure  of  fit  of  the  experimental 

data  to  these  models,  the  mean  deviation  of  the  experimental  points  from  the 

-3  -I 

simple  diffusion  model  above  12  Amagat  is  (all  in  units  of  10  cm  )  +  0.  9 
and  the  root  mean  square  deviation  is  ±3.  For  the  hard  collision  model 
below  12  Amagat  the  mean  deviation  is  +0.  07  and  the  rms  deviation  ±  7.  5. 
Noise  in  the  points  above  12  Amagat  is  primarily  laser  flicker.  Statistical 
noise  is  important  below  12  Amagat. 

The  dashed  line  gives  the  behavior  of  a  Brownian  motion  or  "soft 
collision"  model  (2,4)  for  the  same  region  and  with  the  same  value  of  Dq. 

This  model  assumes  a  small  velocity  change  in  a  single  collision.  If 
DQ  is  chosen  1 0%  lower  for  a  best  fit  to  the  Brownian  motion  model,  the 
mean  deviation  below  12  Amagat  is  -5  and  the  rms  deviation  ±10.  Above  12 
Amagat  the  mean  deviation  is  +3  and  the  rms  deviation  is  ±4.  This  is  a 
noticeably  worse  fit  than  for  the  hard  collision  model,  and  indicates  that 
the  hard  collision  model  is  a  better  approximation  for  hydrogen. 

We  are  presently  extending  these  measurements  to  other  lines  in 
hydrogen  and  deuterium. 

We  would  like  to  acknowledge  very  helpful  discussions  with  Professor 


A.  D.  May. 
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FIGURE  CAPTIONS 


Fig.  1:  Sample  chart  recordings  of  the  Q(l)  R;  man  line  of  hydrogen,  show¬ 
ing  the  dependence  of  linewidth  on  pressure.  The  small  peak  to  the  right  is 
the  Q(0)  line  in  another  order  of  the  interferometer.  The  scanning  speed 
was  0.  025  cm"1 /min  and  the  time  constant  10  sec. 

Fig.  2:  Linewidth  (full  width  at  half  maximum)  of  the  Q(l)  Raman  line  of 
hydrogen  for  forward  and  backward  scattering.  The  Doppler  broadening 
is  destroyed  as  the  density  increases,  and  pressure  broadening  then  be¬ 
comes  dominant.  The  instrumental  linewidth  has  been  subtracted.  The 
theoretical  fits  are  explained  in  the  text. 


FOOTNOTES 


The  line  shape  should  be  Lorentzian,  and  fits  a  Lorentzian  to  within 
the  error  of  the  experiment,  above  12  amagat  in  the  backward  direction 
and  2,  5  Amagat  in  the  forward  direction,  so  that  the  Lorentzian  instrumental 
width  may  be  simply  subtracted.  Below  12  Amagat  the  lineshape  should 
go  smoothly  from  a  Lorentzian  to  a  Gaussian  at  zero  density,  but  our 
noise  at  these  low  densities  is  too  severe  to  get  a  detailed  line  shape.  We 
have  used  a  best  Lorentzian  fit  to  extract  the  linewidth. 


o 
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Frequency  Spectrum  of  Spontaneous  and  stimulated  Line  Narrowing 
Liiect-s  Induced  b\  Laser  Radiation 


M.  S.  Feld  and  A.  Javan 


In  an  earlier  publication,  laser-induced  line  narrow  ingetfect  was 

utilized  in  a  precise  determination  of  isotope  shifts  in  two  Ne  transitions.  A 
<  2) 

recent  paper  also  reports  observation  of  this  effect  and  its  application  in 
studying  some  linewidth  parameters  in  Nt .  These  experiments  study  fluo¬ 
rescence  arising  from  the  low  er  level  of  a  Doppler-broadened  laser  transition 
toamird  level.  [See  Fig.  1(a).]  Viewed  along  the  laser  axis,  the  broad 
fluorescence  line  shape  is  dramatically  influenced  by  the  laser  field.  For  a 
standing-wave  field  detuned  from  its  atomic  center  frequency,  the  laser-induced 
change  signals  appear  as  resonant  increases  in  intensity  over  two  narrow  inter¬ 
vals  symmetrically  located  on  opposite  sides  oi  the  fluorescence  center  frequency. 
The  fluorescence  from  the  upper  laser  level  [  Fig.  1(c)  J  ,  similarly  viewed, 
would  exhibit  narrow  resonant  decreases  in  its  overall  emission  profile. 

The  overall  features  of  this  effect  may  be  described  by  noting  that  the  stand¬ 
ing-wave  laser  field  selectively  interacts  with  atoms  whose  velocities 
Doppier-shift  one  of  its  travelling-wave  components  into  resonance;  this 
produces  changes  in  the  laser  level  populations --an  increase  in  the  lower 
level  population  and  a  decrease  in  the  upper  level  population --over  two 
narrow  intervals  symmetrically  located  about  the  center  of  the  •  elocity 
A  recent  letter^  has  analyzed  the  line -shape  details 


distribution. 
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for  tiir  cascade  case  U  tg.  1< ]  in  terms  oi  two-quantum  transitions 
from  level  2  to  level  1,  are  predirts  differing  widths  for  the  two 
laser- induced  change  signals.  \  similar  line  shape  asymmetry,  de¬ 
scribed  below,  would  appear  in  the  change  signals  from  the  upper  laser  level 
IFig.  1(c)]  .Note  that  in  the  latter  case,  however,  the  0—  1  emission  act  is  an 
inherently  single-quantum  event,  requiring  another  description.  The  sum¬ 
mary  of  a  different  treatment^  based  on  the  density  matrix  formalism  has 
described  the  spontaneous  emission  profile  in  both  cases.  Lineshape  be¬ 
havior  of  similar  origins  has  been  encountered  in  the  microwave  region  in 

experiments  involving  the  interaction  of  two  monochromatic,  classical  fields 

#  -  0^ 

with  a  three-level  system.  ’  ‘  These  considerations  mav  be  readily  extended 
to  the  present  case  by  including  the  Doppler  effect,  and  noting  that  the  spontaneous 
emission  spectrum  from  either  laser  level  lie  vei  Oof  Figs.  1(a)  and  1(c)]  follows 
ihe  emission  line  shape  stimulated  by  a  weak  monchromatic  field  tuned  through 
that  resonance  when  the  population  of  level  1  is  ignored.  In  fact,  recent  experi¬ 
ments  have  studied  the  response  of  a  Doppler-broadened  three-level  svstem  coupl¬ 
ed  to  two  classical  fields.  The  line-shape  theory  describing  these  experiments 
is  directly  applicable  to  important  special  cases  of  the  effect  observed  in 
spontaneous  emission.  This  letter  has  three  intimately  related  objectives:  Th< 
first  is  to  point  out  the  direct  relevence  of  the  detailed  analyses  of  Refs.  7  and 
8  to  this  problem ;  these  also  predict  that  one  of  the  laser-induced  resonances 
will  be  narrower  than  the  other.  Secondly,  we  emphasize  that  experiments 
based  on  stimulated  versions  of  this  effect  demonstrate-  the  different  charac¬ 
teristics  of  the  two  laser- induced  resonances.  Thirdly,  this  letter  gener -iliz.es 
the  theory  in  several  wa  s,  including  importai  power-broadening  effects. 

The  analv  sis  of  Ref.  7  describes  the  third-order  interaction  of  two  mono¬ 
chromatic  fields  with  the  folded  Doppler-broadened  system  of  Fig.  1(b)  in  which 
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l*vis  1  and  2  art  assumed  to  be  closely  spaced.  Note  that  the  analysis  of  a  cascade 

system  {Fig.  1(a)  ]  in  which  the  middle  level  ties  about  halfway  between  the  other  two 

levels  follows  identically.  For  purposes  of  illustration,  let  us  assume  that 

both  transitions  exhibit  gain.  Consider  a  strong  standing-wave  laser  field 

Eg  (Qg)  detuned  by  an  amount  !A9P  +  ^rom  “2,the  Ppilk  2-0  DoPPlcr  re- 

sponse  lFig.  1(b)]  ;here,  7-^=  i  (y .+ Vj).  with  't\  the  natural  decay  rate  of  level  j. 

As  a  weak  monchromatic.  field  (F^ )  is  tuned  through  the  1-0 transition,  it  is  shown 

that  the  broad  Doppler  gain  profile  is  considerable  modified  by  the  presence  of  Eg. 

( 91 

The  resulting  line  shape  is  obtained  from  Eqs.  ( 33d)  and  ( 33e)  ofRef.  7: 


P(«L)  =  G(ot)  f  i  + 


(j  ^ 

w2  * 


eE?  Im{(A2^1)  +  i(721  +  T0)  +  (A9  +  i721  }}  U) 


in  which  Aj  =  atomic  center  frequency  of  .i-0  transition,  and  FL  and 

E°are  frequency  and  amplitude  of  Ei  (FL),  respectively;  G(C^)  is  the  linear 
Doppler  response,  a  slowly  varying  function  ofFl^.  and  £  is  a  constant 
factor  Fo  r  f-  9  above  u9  (1)  predicts  narrow  l.orentzian  responses  of 
widths  Fg  =  7^+70*27q  below  -  ^  and  7.,  above  the  latter  be¬ 

ing  narrower  by  27^  and  independent  of  7q.  In  subsequent  discussion  we 
shall  refer  to  the  "broad"  (F^)  and  "narrow"  (r^-)  resonances,  respectively. 

The  narrow  resonance.  predicted  by  (1)  has  been  observed  and 
fully  verified^) for  Ij.  <  In  these  experiments  levels  1  and  2  are 

tunable  Zeeman  components  of  an  upper  laser  level  connected  to  a  common 
lower  level, 0.  The  monochromatic  fields  are  two  oscillating  laser  modes 
of  fixed  frequencies  determined  by  the  cavity  length.  The  ffert  is  observed 
as  sharp  decreases  m  the  laser  output  as  the  2,  1  level  splitting  is  magnet- 
icallv  tuned.  These  methods,  however,  have  not  been  applied  to  the  ob¬ 
servation  of  rD,  which  requires  absolute  frequence  control  as  in  Lamb-din 

D 


7.  ) 
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experiments  (See  Kef. 


We  report  here  an  extraordinary  manifestation  of  the  broad  resonance, 

3  3 

observed  in  high -resolution  studies  of  the  3p  P^  1  -  3s  Sj  atomic  oxygen 
fine-structure  laser  oscillations  at  #446  A.  The  relevant  fine -structure 
components  of  this  system  consist  of  two  resolved  and  closely  spaced 
Doppler -broadened  transitions  forming  folded  three- level  systems  of  the 
type  shown  m  Fig.  1(b).  For  reasons  unrelated  to  the  present  discussion, 

(i)  the  central  portion  of  each  ol  the  fine-structure  gain  profiles  is  entirely 
depleted  and  appreciable  gain  exists  only  on  the  wings:  and  (ii)  laser  action 
most  readily  occurs  on  the  high-frequency  side  of  the  strongest  (2-0)  fine- 
structure  component.  For  each  laser  mode  oscillating  on  the  2-0  transition, 
two  Lorentzian  holes  of  enormously  different  widths  selectively  deplete  the 
1-0  Doppler  gain  profile:  the  low -frequency  hole,  of  width  2-^+7  (->-y ,,=1 30  MHz,  is 
about  15  times  broader  than  the  high-frequency  hole,  of  width  y^^y0-  9  MHz.  ^  ^ 
Because  of  multimoding  on  the  2-0  transition,  multiple  hole  pairs  are  burnt 
into  the  1-0  transition.  The  broad  low-frequency  holes  overlap,  completely- 
suppressing  laser  action  below  the  1-0  line  center,  and  oscillation  can  occur 
between  the  narrow,  nonoverlapping  holes  above  the  center  frequency.  The 
8446  A  spontaneous  emission  and  laser  output  have  been  studied  photographic¬ 
ally’  with  a  high -resolution  Fabry-Perot  interferometer  with  free  spectral 
range  of  0.  8  cm  *  under  a  wide  range  of  conditions  The  laser  oscillations 
were  obtained  in  argon-oxygen  and  neon-oxygen  mixtures,  using  a  3- 
meter  cavity  with  50  MHz  mode  spacing.  The  operating  pressure  was  kept 
below  1  torr  to  minimize  pressure  broadening  effects.  To  establish  the 
frequency  shifts,  Fabry-Perot  images  of  the  laser  oscillations  and  the 
spontaneous  emission  were  superimposed  upon  the  same  glass  plate  emul¬ 
sion.  In  absence  of  laser  oscillation,  the  spontaneous-emissiomanalysis 
shows,  m  order  of  increasing  frequency”,  a  well-resolved  fine-structure 
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eomponent  (1-0)  with  a  completely  symmetrical  profile,  and  the  strongest 
fine-structure  component  (2-0),  overlapped  on  the  high-frequency  wing  by  the 
third  (weak)  fine-structure  component.  Close  to  threshold,  laser  oscillation 
first  occurs  on  the  high-frequenev  wing  of  the  2-0  component,  where  maxi¬ 
mum  gain  occurs,  as  explained  in  Footnote  11  (ii).  Further  above  threshold 
the  1-0  transition  also  breaks  into  oscillation.  Despite  the  observed 
symmetry  of  the  1-0  profile,  this  concurrent  oscillation  occurs  above  the 
1-0  center  frequency,  an  effect  which  must  follow  from  the  depletion  of  gain 
over  the  entire  low-frequency  wing,  as  described  above.  (It  must  also  be 

pointed  out  that  in  the  observations  reported  here  the  1-0  oscillation  was  at 

(15) 

all  times  close  to  threshold  and  considerably  weaker  than  the  2-0  oscillation.  ) 

Reference  7  treats  the  field  interactions  up  to  third  order.  We  now  outline 
a  different  approach/8^  formulated  in  terms  of  single-  and  double -quantum 
transitions,  which  includes  power-broadening  effects  due  to  a  strong  laser 
field.  The  latter  are  of  considerable  interest:  on  the  theoretical  side,  it  is 
important  to  inspect  whether  or  not  the  strikingly  simple  line-shape  behavior 
is  merely  characteristic  of  a  third-order  culation;  on  the  experimental 
side,  it  is  important  to  know  the  influence  of  power  broadening  on  the  observed 
line  shape.  Consider  the  cascade  system  of  Fig.  1(a).  [The  final  result  will  be 
written  in  a  form  also  valid  for  emission  from  the  upper  laser  level.  Fig.  1(c).  ] 
It  is  desired  to  calculate  the  0-1  emission  spectrum  stimulated  by  the  weak 
travelling-wave  field  E^(D^)  in  the  presence  of  the  strong  standing-wave 
field  E9(I79)  coupled  to  the  2-0  transition.  Specifically  standing-wave  effects 
are  avoided  by  taking  detuned  (  A9  '  >  -^q);  oppositely  propagating 

travelling-wave  components  of  E9  do  not  couple,  and  the  interaction  consists 
of  these  components  independently  coupled  with  Ep 

Consider  an  ensemble  of  atoms  moving  with  given  axial  velocity  v. 


the  f  7  axis  being  defined  by  the  propagation  direction  of  the  travelling- wave 
field  Ej.  In  the  atoms’  rest  frame  the  incident  fields  appear  Doppler- shifted 
to  frequencies  1^(1  -  -)  and  12.,  =  <2.,  (1  -  €  _)  in  which  €  =  +  1  indicates  the 
travelling-wave  component  of  E0  propagating  along  the  +  z  direction,  respec¬ 
tively.  We  now  require  a  solution  in  the  atoms’  rest  frame  in  which  the 
three-level  system  is  coupled  to  one  of  the  travelling-wave  components  of 


the  strong  field  E9  at  frequency  ~  and  to  the  weak  travelling-wave  field 

1 

at  frequency  £2^  The  resonant  interaction  of  two  monochromatic 

fields  with  a  three-level  system  was  treated  in  Ref.  5  for  the  case  of  7.  all 

.1 

equal.  The  perturbation  method  consisted  of  first  obtaining  a  closed  form 
solution  to  the  Schrodingei  equation  for  E°  -  0  and  E^  arbitrary;  and  then 
using  this  result  to  generate  a  solution  valid  to  first  order  in  E”.  When  the 

/t  e  \ 

method  is  extended  to  the  case  of  arbitrary  7.,  1  the  emitted  power  in- 

3 

duced  by  Ej  at  is: 


R  -  iP/ 4fi  I,  and  |i^0  is  the  electric- 

dipole  matrix  element  connecting  levels  i  and  0;  n  is  the  number  of 
background  atoms  with  velocity  v  in  level  3,  i.  e.  the  population  in  absence 
of  the  strong  laser  field.  In  Eq.  (2)  the  n^  coefficient  is  obtained  from 

1  I 

the  2  -*  1  transition  rate  due  to  double-quantum  emission  at  Slj  and 
then  nQ  coefficient,  in  contrast,  is  obtained  from  the  single  -quantum 
emission  rate,  arising  from  0  —  1  transitions  as  modified  by  the  presence 
of  E9,  the  n}  coefficient  results  from  the  reverse  processes,  namely, 
double -quantum  1  —  2  transitions  and  single- quantum  1  -*  0 
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transitions. ^ (See  Ref.  6  )  As  a  check  of  the  detailed  algebra,  Eq.(2) 
lias  also  been  oblained  in  an  independent  calculation  using  tin'  ensemble- 
averaged  density  matrix  to  estimate  tie  induced  polarization  at  an 
approach  equivalent  to  the  one  presented  here.  ' 

It  is  important  to  note  that  (2)  is  entirely  valid  for  the  ease  in  which 

the  yds  are  interpreted  as  decay  rates  arising  front  hard  collisions.  The 

J  ( 5) 

detailed  features  of  the  line  shape  predicted  by  (2)  for  equal  y.  have  been 

J 

fully  verified  in  the  microwave  region  where  Doppler  effect  is  negligible  and 

\ 

the  linewidths  are  entirely  due  to  collision  effects. 

The  complete  emission  spectrum  is  obtained  by  averaging  (2)over  the 

atomic  velocity  distribution  for  €  =  +  1  and  for  t  ~  -  1.  In  the  fully  Doppler- 

broadened  limit  y/D  «  1  (y  ~y. .  and  D=  Doppler  width),  and  for  <*>.  w  , 

1J  l  C» 


Nn  -  Nr 


P(f2  .)  =  G(rt  )  f  1  +  ~  Z  E°2lm I - 1 _ —  +  _ L___  11  ( 3) 

10 1 1  i^ai+ <7^a2) +i-2B  <ai-:j2a2) 

m  which  irB^'io+Si720Q  +  ^t(i':)  ari,1irN  710  +^Y20Q  ‘^(Q  +  1);  Nj  is  the  total 

2 

background  population  of  level  j,  is  a  proportionality  factor  >  0,  Q2  =  i  +  4  Iff  2  I 

o2  y0y2 

and  ^o'^l^l  is  the  power  emitted  at  Slj  by  ihe  Doppler-broadened 

0-1  transition  induced  by  in  the  absence  of  the  laser  field.  Equation  (3)  shows 

the  power  broadening  influence  of  the  laser  field,  which  enters  in  a  remarkably 


simple  way.  Equation  (3)  has  been  written  in  a  form  valid  for  both  cascade 

(a  -  1)  and  folded  {a  =  +  1)  cases,  Figs.  1(a)  and  1(c). 

As  pointed  out  above  and  explained  in  Ref.  16,  the  spontaneous  emission  spectrum 

from  the  upper  or  lower  laser  levels,  as  viewed  along  the  laser  axis,  is  also  given 
by  Eq.  (3)  when  N^is  set  equal  to  0  and  G(£2j)  is  interpreted  as  the  usual  Doppler- 
broadened  spontaneous  emission  spectrum  for  E9  =  0. 

1'ne  discussions  of  Ref.  3  are  consistent  with  weak-field  limit  of  our  treat- 
men:  as  it  applies  to  spontaneous  emission  for  the  special  case  of  Nq=0.  Note, 
for  instance,  that  for  .3  i  «  v  and  in  the  limit  of  y^  -’0,  the  frequency  depen¬ 
dence  of  the  2  -»1  transition  rate,  obtained  from  the  n9coefficient  of  our  Eq.  (2), 


-108- 


would  involve  a  6-function,  becoming  +^2  '  '  “g)  ILgl  ^ the  Ygj3  0  discussion 

of  Ref.  3  is  equivalent  to  averaging  this  distribution  over  velocities. 

To  emphasize  the  significance  of  the  role  played  by  the  background 
atoms  in  level  0,  consider  a  cascade  system  in  which  only  level  0  is  populated 
(i.  e.  Nj  =  N,^=  0).  Then  in  the  rest  frame  of  an  atom,  an  applied  laser  field  at 

i  i 

^2  will  diminish  the  transition  rate  atS2^,  leading  to  two  holes  of  width  r„ 
and  rx  superimposed  upon  the  emission  profile.  [See  Eq.  (3).]  As  dis- 

(5) 

cussed  earlier,  a  0  ~*\  transition  is  an  inherently  single-quantum  event 
and  may  not  be  described  in  terms  of  a  double-  quantum  process  as  in  a  2—1 
transition. 

As  pointed  out  above,  in  the  oxygen^10'  * and  xenon^10^ experiments 
Fg  and  rN  differ  enormously .  In  the  Ne  spontaneous  emission 
experiments  reported  ir  Refs.  1  and  2,  however,  they  are  expected  to  differ 
by  only  about  .30%.  The  observation  of  this  difference  would  require  high 
finesse  Fabry -Perot  analysis  and  good  laser  stability  and  has  not  yet  been 
achieved. 

In  a  veraging  Eq.  (2)  for  the  case  of  finite  72]  a  number  of  cancellations 
occur  in  the  fully  Doppler-broadened  limit  (7/ D^<  1),  leading  to  a  par¬ 
ticularly  simple  expression.  It  is  important  to  point  out  that  such  cancel¬ 
lations  do  not  occur  in  higher  orders  of  7/D.  For  instance,  the  complete 
cancellation  of  in  T^.  which  occurs  in  the  case  of  wj/w2  u  not 

occur  in  the  next  order  of  7/D. 

A  paper  including  complete., algebraic  details  and  additional  discussions  is 
being  submitted  for  publication  .  ; 
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1  \  -tv-w  (H  tin-  intf-n.si?>  <il  1-0  1  •  xer  n.xrili  ••  t i » i r  fu  fct i»>n  of 

ravitv  If  ngfh  would  in-  of  wtft  n  ->t. 

16.  wi  this  cast-  the  time-dependent  a  ave- function  V  is  outuino  1  from  a 

three -level  Schrodinger  equation  to  which  radiative  decav  terms  have 

2  -  i  U  ,t 

oeen  added.  For  ♦  =  e  u  ,  with  u.  the  eigenfunction  of 

level  j  of  energy  *  \V  the  coupled  eouations  arc:  6.  -  Z  (a  -  i  -/  ,  )c 

.  .  .  FHj  Eft)  ‘(Uj-W,):  1  IJ  2  ’J  1  ’ 

m  which  a.  =  -  - e  ,  and  EC)  is  the  sum  of  the  two 

travelling- wave  fields  as  seen  in  the  atoms*  rest  frame 


i  f , 


As  an  example,  for  ar.  atom  in  level  0  sit  initial  tisnc^  ?  -C-Ct-t  t  )  *5 

0  j  0  j  0 

and  the  0  -*  l  transition  rate  at  a  latr  .-ume  t  is  7j  r  (i.  t  )  2.  (St ,  prereeding 

footnote.)  Thus  the  total  stimulated  power  emitted  t> v  background 

1  \  ■> 

atoms  in  1«-  el  0  is  *  n07Q7 ^  !  e,fl,  t  i  '  “  dtQ. 


18.  In  extending  Ivr.  (2)  to  the  spontaneous  emission  rase,  the  ponuiation  of 
level  1,  ,  should  be  set  equal  to  0  and  the  energy  density  of  the  weak 

o  ■>  •  3  i  S  » 

probe  field,  E.  81,  should  be  replaced  by  (hCt  8  v  e  )  d'.T  dS,  where 

1  'll 

frequency  interval  dfi*  «  7  and-  nS  is  1  small  solid  angle  in  the  forw  ard 
direction  {*  7  -  axis)-  Eg.  the  laser  field,  rental  s  in  its  classical  mono¬ 
chromatic  form.  Similar  remarks  apply  to  Eq.  (3);  note,  in  particular, 
that  G  ( )  becomes  the  usual  Doppler-broadened  snec  nurr.  of  the  power 
emitted  spontaneously  into  d  il  dS  with  given  polarization. 


r  IGUR E  C A PT IONS :  The  figure  caption  is  contained  in  the  figure. 
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